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ABSTRACT 


A  combined  experimental  and  theoretical  investigation  of  the 
radiation-induced  impulse  and  stress  generation  characteristics  of 
aluminum  has  been  carried  out  for  cases  where  a  significant 
portion  of  the  blowoff  material  was  in  the  liquid-vapor  state. 
Intense  pulsed  electron  beams  were  used  as  a  source  of  rapid 
(50  nsec)  high-energy-density  loading  to  produce  mixed-phase 
effects  in  aluminum.  Impulse  measurements  were  obtained  with  a 
ballistic  pendulum  for  peak  doses  of  up  to  3800  cal/gm.  Manganin 
gauges  and  the  laser  interferometer  were  used  to  obtain  stress 
time  records  adjacent  to  the  region  of  energy  deposition. 

A  detailed  comparison  of  the  PUFF-V,  GRAY,  and  PHILCO-FORD 
equation-of-state  models  was  carried  out  for  aluminum.  The 
experimental  data  were  compared  with  the  predictions  of  the 
LINEAR  FLUID,  PUFF-V,  GRAY,  and  PHILCO-FC’D  equation-of-state 
models.  The  PHILCO-FORD  equation-of-stat  del  gave  the  best 
data  representation.  The  PUFF-V  model,  which  includes  a  vapor 
state,  and  the  LINEAR-FLUID  model,  which  does  not,  gave  approxi¬ 
mately  the  same  results  for  the  same  input  conditions.  For  the 
experimental  conditions  obtained,  the  impulse  in  aluminum  was 
dominated  by  the  behavior  of  the  equation  of  state  of  aluminum 
in  the  liquid  region. 
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The  laser  interferometer  data  for  the  high-fluence  conditions 
and  the  predictions  of  the  GRAY  equation-of-state  model  suggested 
the  possible  existence  of  a  double  wave  structure  in  the  propagated 
stress  pulse.  This  phenomenon  was  attributed  to  the  recompression 
of  the  liquid  layer  by  the  overlying  vapor  layer  after  the  initial 
drop  to  zero  pressure  in  the  liquid. 
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SECTION  1 
INTRODUCTION 


The  response  of  materials  to  rapid  energy  deposition  is  a 
complex  process  which  depends  on  a  combination  of  physical 
phene nena.  The  rapid  energy  deposition  can  be  the  result  of  a 
nuclear  explosion,  deposition  of  electrons  by  intense  electron 
beams,  or  the  heating  and  compression  of  a  D-T  pellet  by  pulsed 
laser  irradiation  for  controlled  thermonuclear  fusion  applica¬ 
tions.  All  of  these  energy  sources  are  capable  now  or  in  the 
near  future  of  providing  internal  energy  densities  in  excess  of 
the  vaporization  energy  of  most  materials.  The  characteristic 
heating  time  for  these  sources  can  vary  from  a  few  tens  of  pico¬ 
seconds  to  a  few  hundred  nanoseconds. 

The  ability  to  predict  and  understand  the  response  of 
materials  to  rapid  energy  deposition  depends  strongly  on  having 
an  equation  of  state  of  the  material  that  is  able  to  follow  the 
series  of  thermodynamic  states  through  which  the  material  passes 
both  during  and  after  the  rapid  energy  deposition.  The  rapid 
heating  combined  with  high  internal  energy  densities  can  result 
in  pressures  in  the  megabar  range  and  temperatures  on  the  order 
of  up  to  a  few  electron  volts.  When  the  material  expands  adia- 
batically  from  the  initial  high-temperature  high-pressure  state 
after  energy  deposition,  the  density,  pressure,  and  temperature 
drop.  The  equation  of  state  in  this  expanded  volume  region  is 
important  in  determining  the  impulse*  imparted  to  the  material  as 
a  result  of  the  rapid  energy  deposition. 

*As  used  throughout  this  report,  impulse  is  defined  as  change  in 
momentum  per  unit  area. 


Uncertainties  in  the  equation-of-state  models  in  the  expanded 
volume  region  can  lead  to  large  differences  in  impulse  calcula¬ 
tions.  Unfortunately  there  are  no  good  experimental  data  avail¬ 
able  in  the  expanded  volume  region  of  most  materials  on  which  to 
base  an  equation-of-state  model. 

The  objective  of  this  work  was  to  experimentally  and  theo¬ 
retically  investigate  the  impulse  and  stress  generation  charac¬ 
teristics  of  aluminum  that  are  the  result  of  rapid  energy  deposi¬ 
tion  of  internal  energy  densities  in  excess  of  the  vaporization 
enthalpy  of  aluminum  at  one  atmosphere.  For  this  condition 
a  significant  portion  of  the  blowoff  material  was  in  the  liquid- 
vapor  state.  This  work  is  an  extension  of  that  reported  in  DASA- 
2596,  Equation  of  State  Evaluation,  April  1972,  DASA-2475,  Melt 
Dominated  Impulse  Experiments  and  Calculations,  September  1970. 

The  purpose  of  the  program  was  to  investigate  the  contribution  of 
material  in  the  liquid-vapor  state  to  the  material  response  of 
aluminum.  The  experimental  data  were  intended  to  provide  a  data 
base  for  the  evaluation  of  existing  multi-phase  equation-of-state 
models  and  provide  a  basis  for  making  appropriate  modifications 
of  the  multi-phase  equations  of  state  when  indicated  by  the  data. 

Aluminum  was  used  as  a  sample  material  for  several  reasons: 

1)  it  is  a  relatively  easy-to-obtain,  inexpensive  metal  and  2)  it 
has  no  major  solid-solid  phase  transitions,  nor  is  it  particularly 
rate-sensitive . 

Intense  pulsed  electron  beams  were  used  as  the  source  of 
rapid  (50  nsec)  high-energy-density  loading  to  produce  mixed- 
phase  effects  in  aluminum. 


SECTION  2 

EXPERIMENTAL  PROCEDURES  AND  RESULTS 


The  experimental  procedures  and  results  are  presented  in 
this  section.  Description  of  electron  beam  diagnostics  will  not 
be  discussed  in  detail  here  since  an  adequate  account  appears  in 
Reference  1.  However,  a  discussion  of  the  backscatter  correc¬ 
tion  to  measured  total  beam  energy  is  given  in  Appendix  A.  The 
backscatter  correction  was  applied  to  the  data  in  this  program. 

Metallic  beam  guides,  vented  to  reduce  spurious  impulse 
from  anode  debris,  will  be  discussed  here. 

Material  response  data,  consisting  of  impulse,  stress  his¬ 
tories  and  mass  loss,  taken  over  a  peak  dose  range  from  300  to 
3800  cal/gm,  will  comprise  the  bulk  of  this  section. 

2.1  BEAM  TRANSPORT  AND  SPURIOUS  IMPULSE 

High-current  electron  beams  generated  by  Physics 
International's  Model  738  Fulserad  were  used  to  irradiate 
aluminum  samples.  Dynamic  loads  were  induced  in  the  test 
material  by  rapid  deposition  of  electrons  (50  nsec  FWHM  of  beam 
power  curve)  in  the  front  surface  of  targets.  Peak  energy  dose 
levels  were  varied  from  300  to  3800  cal/gm  by  changing  the 
intensity  (fluence)  of  200  to  400  keV  mean  electron  energy 
beams.  Active  beam  diagnostics,  described  in  Reference  1,  were 


used  to  characterize  the  electron  beam  environment  and  to 
compute  energy  deposition  as  a  function  of  material  depth. 

Peak  dose  in  the  material  was  varied  by  changing  the  elec¬ 
tron  beam  fluence.  The  method  of  varying  fluences  for  a  speci¬ 
fied  mean  electron  energy  (and  hence,  for  fixed  total  beam 
calories)  consists  of  concentrating  the  total  calories  in  the 
beam,  H,  into  different  uniform  irradiation  areas,  A.  The  test 
matrix  specified  the  mean  electron  energy  and  the  fluence.  For 
a  given  mean  electron  energy  the  total  machine  output  in 
calories  is  fixed.  Therefore,  the  required  irradiation  area  was 
then  determined  by  A  =  H/?. 

Metallic  beam  guides  and  gas  focusing  were  used  for  beam 
transport  and  fluence  control.  Background  gas  pressure  in  the 
electron  beam  drift  chamber  governs  the  magnitude  of  the  net 
current  and  hence  the  degree  of  beam  pinching  due  to  self- 
magnetic  fields.  Metallic  guide  cones  were  used  to  shape  the 
gaussian  fluence  at  the  anode  into  a  uniform  distribution  over 
the  sample  placed  at  the  cone  exit.  Two  beam  guide  geometries 
were  used  to  achieve  a  broad  range  of  fluences.  The  beam  guide 
for  low-to-medium  fluences  (approximately  10  to  50  cal/cm2)  was 
a  20-cm  long,  straight,  stainless-steel  cone  with  slotted  walls 
(Figure  1).  The  high-fluence  beam  guide  (Figure  2)  consisted 
of  a  short,  solid,  graphite-lined,  stainless-steel  cone  with  a 
curved,  slotted  pipe  attached  to  its  exit. 


The  slots  in  the  beam  guide  walls  were  designed  to  minimize 
spurious  impulse  from  anode  debris,  the  source  of  which  will  be 
discussed  below.  The  location  of  slots  in  the  low-to-medium 
fluence  beam  guide  (Figure  1)  needs  no  explanation.  The  func¬ 
tion  of  the  slotted  pipe  attached  to  the  exit  of  the  high- 
fluence  beam  guide  was  to  bend  the  beam  (through  approximately 
45  degrees)  away  from  its  original  axis  while  allowing  anode 
debris,  collimated  by  the  cone,  to  escape  through  the  slots  in 
the  pipe  walls. 

The  following  experiments  were  performed  to  demonstrate  the 
negligible  anode-debris  contribution  to  impulse. 

The  anode  of  the  738  Pulserad  is  a  1/4-mil  sheet  of 
aluminized  Mylar.  Anode  material  covering  the  injection  area  of 
the  guide  cones  is  vaporized  on  every  shot  because  of  energy 
deposition  in  the  aluminized  Mylar.  Half  of  the  vaporized  mass 
of  anode  material  expands  into  the  drift  chamber  and  hence  into 
the  guide  cone.  If  the  0.01  gram  of  anode  debris,  moving  at 
sufficient  velocity,  were  to  impact  the  sample  in  a  pendulum  bob 
at  the  end  of  the  guide  cone,  the  pendulum  could  record  spurious 
impulse.  Venteci  (i.e. ,  slotted)  beam  guides  have  been  used  to 
demonstrate  that  the  anode-debris-induced  impulse  level  is  less 
than  100  taps. 

For  the  low-to-medium  fluence  beam  guide,  the  beam  environ¬ 
ment  chosen  was  one  for  which  quartz  gauge  data  for  aluminum  was 
available.  These  traces  typically  show  triangular  stress  pulses 
with  less  than  0.2-psec  duration  (base-to-base) .  Guided  by 
these  data,  a  70-mil-thick  aluminum  disk  (flier)  was  attached 
with  a  thin  film  of  vacuum  grease  to  the  rear  surface  of  a 
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larger  diameter,  50-mil- thick  aluminum  disk  (sample).  This 
"sandwich"  was  placed  in  an  annular  pendulum  bob,  so  that  only 
the  sample  was  held,  while  the  flier  was  free  to  fly  off 
unimpeded.  The  array  was  placed  at  the  end  of  the  straight, 
slotted  guide  cone  vised  in  low-to-medium  fluence  experiments 
(Figure  1).  Each  time  beams  were  fired  into  such  arrays,  the 
fliers  separated,  while  the  bobs  holding  the  samples  showed  zero 
deflection. 

The  separation  time  of  the  flier  subsequent  to  electron 
beam  deposition  in  the  sample  (sample  thickness  plus  twice  the 
flier  thickness)  was  chosen  to  be  1.08  Msec.  At  the  time  of 
separation,  the  flier  became  a  momentum  trap  capable  of  con¬ 
taining  more  than  the  measured  stress  pulse  under  similar  beam 
conditions.  Thus,  1.08  nsec  after  deposition  time,  the  sample 
at  the  cone  exit  was  in  position  to  receive  any  spurious  impulse 
and  record  it  as  a  pendulum  deflection.  No  deflections  were 
recorded. 

Arrival  of  anode  debris  at  the  sample  in  less  than  1.08  nsec 
would  require  40  times  the  available  energy  in  the  electron  beam 
and  would  result  in  impulse  levels  an  order  of  magnitude  greater 
than  the  highest  ever  recorded  at  Physics  International.  There¬ 
fore,  we  concluded  that  no  spurious  impulse  due  to  anode  debris 
can  occur  less  than  1.08  nsec  after  beam  deposition  time.  Also, 
no  measurable  spurious  impulse  occurred  after  1.08  nsec,  at  the 
exit  of  the  20-cm,  slotted  guide  cone  used  for  low-to-medium 
fluence  material  response  experiments. 
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The  absence  of  spurious  impulse  due  to  anode  debris  could 
not  be  demonstrated  for  the  high-fluence  beam  guide  in  the  same 
way  as  for  the  low-to-medium  fluence  guide  cone.  Since  there  is 
reason  to  expect  long  time  impulse  generation  in  aluminum  for 
the  liquid-vapor,  mixed-phase  regime,  no  criteria  existed  for  a 
momentum  trap  thickness  such  that  all  the  real  impulse  would  be 
carried  off  by  the  flier  well  before  possible  anode  debris 
arrival.  The  alternative  demonstration  consisted  of  firing  a 
beam  (less  than  1000  cal/gm  peak  dose)  into  a  graphite  sample  in 
the  pendulum  at  the  exit  of  the  beam  guide.  No  pendulum  deflec¬ 
tion  was  observed,  suggesting  that  spurious  impulse  due  to  anode 
debris  through  the  high-fluence  beam  guide  is  less  than 
measurable  by  the  ballistic  pendulum  used  in  this  program. 

2 . 2  IMPULSE  MEASUREMENTS 


The  impulse  data  presented  here  were  obtained  using  the 
ballistic  pendulum  technique  described  in  Reference  1. 


The  linear  fluid  model  of  impulse  generation  in  aluminum 
was  shown  to  be  valid  in  the  melt-dominated  regime  (Reference  2) . 
It  states  that 

I  =  0.04186  ~  <j>  exp 
"*k 


where 


I  =  impulse  in  ktap 

=  bulk  sound  speed  (C^  =  >/k/p0)  in  cm/usec 

r  =  Gruneisen  coefficient 

2 

4>  =  fluence  in  cal/cm 
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t  =  deposition  time  in  psec 

R  =  electron  range  in  cm 

X  =  adiabatic  bulk  modulus  in  mbar 

3 

PQ  =  density  in  gm/cm 

Accordingly,  the  generated  impulse  depends,  not  only  on 
fluence,  but  on  deposition  time  and  electron  range  as  well,  at 
least  in  the  melt-dominated  regime.  The  body  of  impulse  data, 
in  this  program,  was  collected  over  a  200  to  400  keV  mean  elec¬ 
tron  energy  spread  with  a  corresponding  spread  of  electron 
ranges  in  the  test  samples.  Furthermore,  since  fluence  control 
was  achieved  with  gas  focusing,  the  variation  of  background  gas 
pressure  introduced  varied  degrees  of  beam- front  erosion 
(Reference  3).  Consequently,  the  beam  pulse  (deposition  time) 
at  the  sample  location  varied  according  to  the  gas  pressure  used 
to  achieve  the  fluence.  In  short,  each  measurement  of  impulse 
can  be  expected  to  have  its  individual  fluence,  deposition  timer 
and  electron  range.  In  these  circumstances,  the  linear  fluid 
model  for  melt-dominated  impulse  generation  predicts  an  artifi¬ 
cial  scatter  of  the  data  when  impulse  is  plotted  versus  only 
fluence — without  regard  to  dependences  on  electron  range  (mean 
electron  energy)  and  deposition  time. 

The  preceding  discussion  was  intended  to  support  the  view 
that  a  comparison  of  model  predictions  with  impulse  data, 
generated  in  electron  beams,  using  impulse  versus  fluence  plots 
as  the  arena  of  comparison  is  not  appropriate.  Such  compari¬ 
sons  must  be  done  on  a  point-by-point  basis  since  each  impulse 
measurement  can  have  its  individual  set  of  loading  conditions 
(e.g.,  fluence,  energy  deposition  range  and  deposition  time)  on 
which  impulse  generation  may  depend.  Therefore,  a  plot  pf 
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impulse  versus  fluence  for  the  electron  beam  data  does  not  give 
a  valid  representation  of  the  scatter  in  the  data,  and  such 
gr  .phs  can  only  be  used  to  present  general  trends. 

The  data  summary  is  shown  in  Table  1,  and  the  notations 
employed  in  it  are  as  follows : 


<E>  =  Mean  electron  energy  (keV)  calculated  as 
described  in  Reference  1. 

<0>  =  Mean  angle  of  incidence  (degrees)  of  the  beam 
electrons  from  the  calculation  in  Reference  3. 

x  =  Full  width  at  half  maximum  (FWHM)  of  the  power 
curve . 

H  =  Total  beam  calories  impinging  on  target. 

—  2 

A  =  Crater  (mass  removal)  area  in  the  target  (cm  ) . 

_  2  _ 

4>  =  Average  fluence  (cal/cm  )  defined  as  H/A. 

A (MV)  =  Momentum  change  (kilo- dyne -sec)  imparted  to  the 
ballistic  pendulum  in  which  the  target  is 
mounted . 

I  =  Impulse  per  unit  area  (ktap)  defined  as  A(MV)/A. 

Am  =  Mass  (gm)  lost  by  the  target. 


Parentheses  indicate  approximate  values.  The  peak  dose  in 
the  target  was  estimated  by  the  product  of  fluence,  (<f>)  ,  and  the 
normalized  peak  dose  characteristic  of  <E>  and  <0>.  It  was 
demonstrated  in  Reference  4  that  <E>  and  <0>  are  sufficient  to 
identify  the  appropriate  normalized  energy  deposition  profile  in 
a  material. 
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2 . 3  STRESS  MEASUREMENTS 


The  laser  interferometer  and  the  manganin  gauge  techniques 
which  were  used  to  record  stress  histories  at  the  rear  surface 
of  samples  are  described  in  Appendices  B  and  C,  respectively. 

The  geometry  of  the  manganin  foil  gauges  is  illustrated  in 
Figure  3. 
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Figure  3  Manganin  gauge  geometry. 

The  gauge  package  was  placed  inside  the  bob  of  a  ballistic 
pendulum  in  order  to  make  simultaneous  measurements  of  impulse 
and  stress  (see  Appendix  C) .  Five  data  points  were  taken  at 
peak  doses  of  approximately  2000  cal/gm.  A  typical  gauge  signal 
is  shown  in  Figure  4.  Qualitatively,  the  stress  histories  con¬ 
sisted  of  sharp,  triangular  spikes  of  approximately  200  to 
300  nsec  followed  by  a  long,  shallow  tail  which  lasted  for  the 
entire  gauge  readtime  (10  psec) . 
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The  millivolt  level  output  signal  of  the  manganin  gauge 
required  extensive  shielding  from  the  rf  environment  in  the 
electron-beam  drift  chamber.  This  required  a  sample  thickness 
of  0.125  inch.  Since  electron  ranges  for  the  irradiations  per¬ 
formed  in  this  program  are  typically  0.020  inch,  the  stress 
measurements  were  made  well  outside  the  deposition  region  and 
the  stress  pulse  was  attenuated. 

Quantitative  interpretation  of  manganin-gauge-measured 
stress  histories  were  further  complicated  by  a  somewhat  non- 
reproducible  hysteresis  (approximately  15  percent)  and  two- 
dimensional  effects  resulting  in  gauge  stretching.  Both  flaws 
are  in  the  direction  of  higher  output  from  the  gauge. 

A  comparison  of  measured  impulse  with  the  integral  under 
the  stress  curve  showed 

Measured  impulse  <  J ’ o(t)  dt. 

The  laser  interferometer  (see  Appendix  B)  was  used  to 
measure  rear-surface  velocity  histories.  The  rear-surface 
velocity  histories  were  then  converted  to  stress  incident  on  the 
rear  surface  by  using  the  Hugoniot  of  6061  alum'.num.  The 
displacement  mode,  applicable  for  low-amplitude  signals,  was 
used  on  thick  (0.5  inch)  samples.  The  velocity  mode,  better 
suited  for  high-  amplitude  pulses,  was  used  with  thin  (0.32  inch) 
samples.  In  both  cases,  the  laser  beam  was  reflected  from  the 
polished  rear  surface  of  the  samples.  Before  the  signals  from 
the  rear  surface  were  terminated  by  rear-surface  spall  (loss 
of  reflecting  surface),  high-resolution  stress  histories  of  the 
thermomechanical  spike  were  recorded  with  the  velocity  mode. 
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The  laser  interferometer  data  is  shown  in  Figures  5  through  11. 
The  elastic  precursor  (at  approximately  7  kbar)  is  clearly  seen 
in  Figures  5  and  6.  The  long-duration,  low-amplitude  tails  were 
recorded  on  separate  shots  at  the  rear  surface  of  thick  samples. 
The  thick  sample  stress  history  in  Figure  8  shows  the  arrival  of 
the  thermomechanical  spikes  after  reflection  from  the  front 
surface.  The  thick  sample  stress  histories  are  undoubtedly  no 
longer  one-dimensional  even  though  their  included  areas  (impulse) 
are  conserved  in  the  two-dimensional  distortion.  Figure  11  shows 
a  composite  of  the  thermomechanical  (thin  sample)  stress  to¬ 
gether  with  the  long-duration  (thick  sample)  stress.  This 
composite  has  qualitative  similarity  with  stress  histories 
obtained  from  manganin  gauge  measurements,  but  exhibits  signifi¬ 
cantly  higher  resolution. 


In  summary,  manganin  gauges  have  been  used  to  measure 
sample  in-material  stress  over  10  ysec  duration  while  simul¬ 
taneously  measuring  impulse  by  positioning  the  gauge  in  a 
pendulum  bob.  This  stress  measurement  technique  yielded  only 
qualitative  information.  Quantitative  interpretation  was 
hindered  by 

1.  Gauge  hysteresis  (non-reproducible) 

2.  Gauge  stretching  (two-dimensional  effects) 

3.  Long  rise time  (equilibration  in  the  gauge  lamina) 

The  laser  interferometer  has  proven  to  be  a  better  tech¬ 
nique  for  recording  stress  histories  since  it  lacks  the 
objectionable  features  inherent  in  manganin  gauge  measurements. 
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Figure  5  Stress  historv  obtained  with  laser  interferometer 
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Figure  o  Stress  historv  obtained  with  laser  interferometer. 
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2.4  MASS  REMOVAL  DATA 

The  existing  aluminum  data  point  out  a  variety  of  crater 
appearances.  A  quick  scan  of  the  data  indicates  that  craters 
vary  in  appearance  with  increasing  fluence  for  a  fixed  mean 
electron  energy.  The  observed  variation  is  as  follows: 

1.  Low  fluence:  craters  are  shallow  with  shiny , 
granular  textured  floor. 

2.  Medium  fluence:  craters  are  deeper  with  dull 
gray,  smoother  textured  floor,  occasionally  dotted  by 
a  few  small  droplets. 

3.  High  fluence:  craters  are  obscured  by  some  liquid 
"splash"  (radially  outward) . 

4.  Highest  fluence:  craters  are  heavily  obscured  by 
"splash,"  but  occasionally  show  deep  crater  floors 
with  large  droplets. 

The  four  descriptions  refer  to  irradiation  at  approximately 
200  to  300  keV.  High- fluence  irradiations  at  about  600  keV 
result  in  deep,  clean  craters  with  shiny,  granular  textured 
floors.  Descriptions  1,  3,  and  4  are  illustrated  in  Figure  D-l. 

In  an  effort  to  account  for  the  diversity  of  crater 
appearances,  the  existing  data  which  cover  a  wide  fluence  (and 
dose)  range  were  analyzed  according  to  crater  appearance,  depth 
of  material  removal,  melt  depth,  electron  range,  and  peak  dose. 
This  analysis  revealed  definite  trends  which  led  to  a  plausible 
explanation  of  crater  appearances  in  terms  of  material  removal 
mechanisms . 
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The  above  analysis  required  an  additional  category  for  the 
data.  It  became  necessary  to  introduce  groupings  according  to 
whether  or  not  the  propagated  stress  was  allowed  to  reflect  from 
the  sample  rear  surface.  Sample  thicknesses  varied  over  the 
data  set  from  0.08  to  0.50  inch  and  some  of  the  thin  samples 
had  momentum  traps  attached  to  their  rear  surfaces .  The  cate¬ 
gory  of  samples  with  no  reflected  stress  consisted  of  thin 
samples  for  which  momentum  traps  "carried  away"  the  propagated 
stress  pulse. 

Based  on  the  criteria  described  above,  all  the  data  were 
sorted  into  two  categories,  each  of  which  consisted  of  two  sub¬ 
groups  according  to  whether  or  not  rear  surface  stress  reflec¬ 
tion  occurred.  Category  A  labeled  samples  with  clear  craters — 
descriptions  1  and  2  above;  category  B  consisted  of  samples 
whose  craters  were  obscured  by  splash — descriptions  3  and  4 
above. 

Before  discussing  trends  established  by  the  data,  it  is 
helpful  to  consider  some  simple  models  of  front  surface  mass 
removal — in  particular,  front  surface  spall  in  a  semi-infinite, 
perfectly  elastic  slab  of  homogeneous  metal.  According  to  the 
derivations  in  Appendix  D: 

1.  The  shortest  distance  from  the  front  surface  at 
which  tension  reaches  its  maximum  value  is  exactly 
half  of  the  electron  range  in  the  material. 

2.  For  finite  spall  strengths,  spall  will  proceed  in 
planes  of  equal  thickness  so  that  the  maximum  spall 
depth  is  less  than  or  (at  most)  equal  to  half  of  the 
electron  range. 


/ 
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Consequently,  when  the  idealized  solid  is  irradiated  with 
an  energy  density  below  melt,  the  material  removal  depth  due  to 
front  surface  spall  will  be  less  than  or  equal  to  half  the 
electron  range,  i.e., 

dR  =  removal  depth  <  half  electron  range  =  R/2 

The  exact  material  removal  depth  will  depend,  of  course,  on  the 
fluence.  The  appearance  of  the  crater  will  be  typical  of  solid 
spall,  clean  and  with  a  granular  texture. 

When  the  irradiation  dose  is  increased  above  melt  but  below 
vaporization  so  that  the  solid  is  melted  to  exactly  half- 
electron  range,  the  same  rule  prevails.  The  exception  comes  in 
the  form  of  thinner  spall  layers  which  allow  the  removal  depth 
to  approach  half  range  more  closely,  i.e.,  dr  »  R/2. 

In  this  case,  all  the  liquid  will  be  removed  to  the  melt 
depth.  The  crater  should  still  look  essentially  clean  but  the 
floor  should  no  longer  resemble  solid  spall.  Its  appearance 
could  be  considerably  smoother  due  to  solidification  of  the 
remaining  material  at  the  incipient  melt  dose. 

The  intermediate  case  occurs  when  irradiation  is  such  that 
the  melt  depth  is  less  than  half  range.  The  liquid  proceeds  to 
spall  to  the  melt  depth.  If  the  tension  produced  by  reflection 
from  the  new  free  surface  at  the  melt  depth  (see  Appendix  D) 
exceeds  the  solid  spall  strength,  spall  will  continue  into  the 
solid  material  but  will  still  not  exceed  half  range,  i.e., 
dR  <  R/2.  The  crater  will  again  be  clean  of  any  liquid  splash 
and  will  have  texture  typical  of  solid  spall. 
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The  crater  appearance  and  irradiation  conditions  discussed 
analytically  above  match  the  data  in  Category  A  for  the  case  of 
no  significant  rear  surface  reflection  of  the  propagated  stress 
(equivalent  to  a  semi-infinite  slab).  Figure  12  shows  a  com¬ 
parison  of  half  range  with  measured  removal  depth  for  that 
data.  It  is  evident  that,  within  experimental  error , 

Removal  depth  <  half  range 

for  all  the  data  in  this  category.  Figure  13  demonstrates  that 
for  some  of  that  data,  melt  depth,  removal  depth,  and,  hence, 
half  range  were  equal. 

If  the  propagated  stress  is  allowed  to  reflect  in  tension 
from  the  rear  (free)  surface  of  a  thin  sample  and  if  the  tension 
exceeds  the  spall  strength  of  the  material  at  that  surface, 
after  rear-surface  spall  and  propagation  back  to  the  front 
surface,  additional  mass  will  be  removed.  The  effect  of  such 
(secondary)  spall  will  be  to  increase  mass  loss  and  removal 
depth  but  will  not  significantly  increase  momentum  imparted  to 
the  sample.  The  particle  velocities  in  the  secondary  spall,  as 
well  as  in  the  remaining  sample,  will  be  in  the  same  direction 
and  the  effect  of  secondary  spall  will  be  to  separate  the  two 
masses,  "trapping"  lower  particle  velocities  in  the  spall. 

Data  in  Category  B,  because  of  its  very  nature,  did  not 
lend  itself  to  accurate  measurements  of  removal  depths.  Liquid 
splashes  of  various  extent  obscured  the  crater.  Furthermore, 
peak  doses  for  much  of  the  data  in  this  category  exceeded 
incipient  vaporization  energy  and  the  partially  vaporized 
material  can  no  longer  be  considered  to  be  linearly  elastic.  It 
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Figure  12  Half- range  versus  measured  removal  depth. 


Figure  13  Melt-depth  versus  measured  removal  depth. 


can  be  reasoned  that  the  pressure  resulting  in  vapor  at  the 
front  surface  can  overcome  the  tension  buildup  in  the  remaining 
material,  in  which  case  the  simple  model  used  above  no  longer 
applies.  However,  that  model  can  still  be  of  use  in  under¬ 
standing  how  splash  occurs  for  irradiation  below  incipient 
vaporization  doses. 

If,  for  example,  the  melt  depth  exceeds  half  range,  front- 
surface  spall  will  remove  liquid  down  to  half  range  leaving  a 
liquid  layer  behind.  If  secondary  spall  is  excluded,  that 
liquid  layer  will  stay  on.  Two-dimensional  radial  relief  can 
result  in  radial  components  of  tension  resulting  in  the  migra¬ 
tion  of  the  remaining  liquid  layer  away  from  the  center  of  the 
crater.  If  solidification  of  the  liquid  occurs  during  this 
process,  the  "splash"  appearance  observed  experimentally  might 
result. 

If  secondary  spall  is  allowed  to  occur  while  the  liquid 
layer  has  begun  localized  solidification,  the  final  appearance 
of  the  crater  may  show  resolidified  liquid  instead  of  the  radial 
splash.  This  condition  will  be  accompanied  by  removal  depths 
far  in  excess  of  half-range.  Examples  of  this  crater  descrip¬ 
tion  are  found  among  the  data. 

As  described  in  the  analysis  section  of  this  report,  the 
experimental  laser  interferometer  data  and  the  predictions  of 
the  advanced  multi-phase  equation-of-state  models  both  suggest 
that  not  all  the  impulse  was  delivered  promptly  (~  100  nsec)  for 
the  high  fluence  and  highest  fluence  experiments.  This  has  been 
attributed  to  the  recompression  of  the  liquid  layer  by  the  over- 
lying  vapor  and  vapor- liquid  layers  after  the  initial  drop  to 
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zero  pressure  in  the  liquid  layer.  This  may  also  account  for 
the  final  appearance  of  the  craters  which  exhibit  the  "splash" 
characteristics  at  the  higher  fluence  levels. 

In  summary,  the  observed  variety  of  crater  appearance  is 
consistent  with  the  front-surface  spall  mechanism,  the  solid- 
to- liquid  phase  transitions,  and  the  predictions  of  the  advanced 
equation-of-state  models.  Data  which  lend  themselves  to 
quantitative  analysis  confirm  this.  When  significant  vapor 
generation  occurs  at  the  front  surface  the  pressure  due  to  vapor 
would  tend  to  suppress  front-surface  spall  of  the  liquid.  If 
radial  relief  or  the  recompression  of  the  melted  liquid  layer 
takes  place,  the  final  crater  appearance  will  exhibit  the 
"splash"  characteristic  commonly  observed  at  higher  doses. 

2.5  SUMMARY  OF  EXPERIMENTAL  PROCEDURES  AND  RESULTS 

Intense  relativistic  electron  beams  were  used  to  achieve 
rapid,  in-depth  heating  in  aluminum  samples.  Impulse  and  stress 
data  were  collected  over  a  peak  dose  range  of  300  to  3800  cal/gm. 

Simultaneous  impulse  and  in-material  stress  measurements 
were  obtained  using  manganin  piezoresistive  transducers  in  con¬ 
junction  with  a  ballistic  pendulum.  Preliminary  stress  data 
were  obtained  using  a  laser  interferometer  in  the  rear-surface 
velocity  and  displacement  modes.  The  laser  technique  demon¬ 
strated  high  signal  resolution. 
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Mass  removal  mechanisms  have  been  examined  and  compared 
with  crater  depths  and  appearances.  Quantitative  agreement  was 
found  in  the  low-to-medium  dose  range  (below  vaporization 
energy) .  Qualitative  interpretation  of  crater  appearances  was 
given  for  the  effects  observed  in  the  higher  dose  range  (above 
vaporization  energy) . 
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SECTION  3 

ANALYSIS  OF  MIXED-PHASE  IMPULSE  GENERATION 

3.1  INTRODUCTION 

The  previous  sections  have  described  the  experimental 
techniques  and  difficulties  encountered  in  obtaining  material 
response  data  on  aluminum  in  the  mixed  phase  regime  (liquid- 
vapor)  .  In  the  program  that  preceded  this  effort  (Reference  2) 
emphasis  was  placed  on  phenomena  found  in  the  melt-dominated 
regime . 

The  objective  of  the  analytical  portion  of  this  program 
was  to  compare  the  experimental  data  with  predictions  using 
existing  multi-phase  equations  of  state.  The  investigation 
was  to  determine  the  importance  of  phase  changes  and  the  input 
parameters  of  the  models  on  various  material  response  character¬ 
istics  ,  such  as  stress  generated  and  impulse  imparted  to  the 
sample.  The  analytical  work  was  also  to  determine  whether  a 
complex  mixed-phase  equation-of-state  model  was  needed  to 
predict  the  observed  material  response  or  whether  a  simpler 
but  more  physical  model,  such  as  suggested  in  Reference  2, 
would  do  equally  well  over  a  specified  range. 

One  of  the  main  results  of  the  preceding  investigation 
(Reference  2)  on  melt-dominated  impulse  was  that  both  the 
experimental  impulse  and  stress  histories  were  predicted  by  a 
simplified  model  which  only  considered  the  solid  and  liquid 
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phases.  A  simplified  formula,  given  in  Equation  1  was  found 
to  predict  the  impulse  for  a  number  of  metals  in  the  melt- 
dominated  regime.  This  simplified  expression  was  found  to  be 
valid  for  energy  densities  well  above  incipient  vaporization  at 
atmospheric  pressure.  A  theoretical  upper  limit  was  established; 
however,  an  experimental  upper  limit  on  the  maximum  energy 
densities  for  which  the  simplified  expression  was  valid  was  not 
established.  This  program  attempts  to  establish  this  upper  limit 
on  energy  densities  for  which  the  simplified  linear  fluid  model 
is  applicable. 

The  experimental  portion  of  this  program  has  obtained  peak 
doses  up  to  approximately  3800  cal/gm  which  is  well  above  the 
energy  density  (3350  cal/gm)  required  to  fully  vaporize  aluminum 
at  atmospheric  pressure.  Because  of  the  high  peak  dose  obtained 
experimentally  in  this  program,  it  was  deemed  advisable  to 
compare  the  data  with  the  predictions  of  the  latest  multi-phase 
equation-of-state  models. 

This  section  contains  a  brief  description  of  the  two  multi¬ 
phase  equation-of-state  models ,  the  PHILCO-FORD  and  GRAY 
equation-of-state  models.  The  experimental  results  are  compared 
with  the  predictions  of  the  following  equation-of-state  models 
in  Section  5; 

a.  PUFF 

b.  LINEAR  FLUID 

c.  GRAY 

d.  PHILCO-FORD 
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3.2  PHILCO-FORD  EOS 

The  PHILCO-FORD  (P-F)  equation-of-state  model  (Reference  5) 
was  developed  by  the  Aeronutronic  Division  of  Philco-Ford 
Corporation  for  the  DNA-sponsored  PREDIX  program.  This  model 
was  developed  for  materials,  such  as  metals,  which  melt  rather 
than  sublime  and  whose  vapor  is  composed  of  atomic  species.  The 
model  includes  provisions  for  treating  the  liquid-solid  and 
liquid-vapor  mixed-phase  regions  as  well  as  the  single-phase 
solid,  liquid,  and  vapor  regions. 

What  follows  is  a  brief  description  of  how  the  P-F  equation- 
of-state  model  treats  each  phase  region  and  constructs  the  various 
phase  transition  lines . 

3.2.1  Vapor-Phase .  In  the  vapor  phase  a  generalization  of 
the  van  der  Waals  equation  is  used.  The  formalism  that  is 
followed  was  first  given  by  Hirschfelder ,  et  al.,  (References  6 
through  8)  for  describing  the  vapor  phase  of  the  noble  gases. 

The  expression  used  is  given  by 

(P  +  A(T)/V2  +  A'(T)/V3)(V  -  B+B'/V)  =  RT  (2) 

where  A  and  B  in  van  der  Waals  equation  are  now 

A  =  A  (T)  +  A'  (T)/V  (3a) 

B  =  B  -  B'/V  (3b) 

with  B  and  B'  being  constants. 
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The  expression  is  then  non-dimensionalized  by  introduction 

of  the  reduced  variables,  p  =  P/P  ,  v  =  V/V  ,  and  t  =  T/T  , 

c  c  c 

where  the  subscript  c  refers  to  values  of  the  thermodynamic 
variables  at  the  critical  point.  With  the  introduction  of  the 
critical  compressibility  ratio,  Zc  =  PcVc/R  T^,  (Equation  2) 
reduces  to 


p  =  pt/Zc  (l  -  bp  +  b'p2)  -  a  (t)  p2  -  a*  (t.)p3  (4) 


where  a(t)  =  A(T)/PcVc2,  a'  (t)  =  A*  (T)/PcVc,  b  =  B/Vc,  b'  =  B*/^2, 

and  p  =  1/v.  The  relationships  between  the  various  parameters 

are  obtained  by  the  requirement  that  the  critical  isotherm  has  a 

point  of  inflection  at  the  critical  point.  In  terms  of  Equation 

2  2~ 

4  this  means  that  (3p/3p)t  and  3  p/3p  ^ are  zero  at  the  critical 
point,  p  =  p  =  t  =  1.  Applying  this  with  the  assumption  that 
a'  (1)  =  0  gives 


a(l)  =  8 


(5a) 


b  =  (382  -  66  -  l)/3 (36  -  U 


(5b) 


b'  =  (8  -  3)/  (36  ~  1) 


(5c) 


where  8  is  related  to  the  critical  compressibility  ratio  by 

Zc  =  6(36  -  l)/(6  +  l)3  (5d) 

Thus  b  and  b'  are  known  once  8  is  determined.  The  determination 
of  8  is  described  in  the  section  that  discusses  the  vapor-liquid 
mixed-phase  region. 

The  following  temperature-dependent  forms  for  a(t)  and  a'(t) 
are  suggested. 
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a(t)  =  kQ  +  ($  -  kQj/t  =  kQ  +  k^t  (6a) 

a'(t)  =  \  (l  -  kQ  +  26  -  a) (t-l/t )  =  k2(t-l/t)  (6b) 

The  parameter  kQ  is  an  input  parameter  and  a  is  determined 
from  considerations  of  the  vapor- liquid  mixed  phase  region  by 
noting  that  a  =  Op/3t)p,  evaluated  at  the  critical  point. 

The  internal  energy  of  the  vapor  phase  is  obtained  by  using 
the  thermodynamic  relationship : 


and  integrating  the  equation  from  zero  density  along  an  isotherm 

P 

E(v.t)  -  E0(T)  =  -Zc  RT^  ^  |t(|E.)0  -  P  K 

-  2c  OTc  |-(ko  +  +  V2/4]  <8) 

where  EQ(T)  is  the  internal  energy  at  the  temperature  T  and  zero 
density.  The  form  for  E(T)  that  is  chosen  is  that  of  an  ideal 
monatomic  gas  with  the  contribution  to  the  internal  energy  from 
electronic  excitations  assumed  to  be  zero.  With  these  assumptions 
Eq(T)  becomes 

EQ(T)  =  £q  +  CVT  (9) 

where  Eq  is  the  reference  energy  at  0°K. 
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The  equation-of-state  for  the  vapor  phase  is  then  given  by 
Equations  4,  8,  and  9.  Given  the  input  Variables  V  and  E,  the 
system  of  equations  is  solved  in  order  to  obtain  the  pressure. 
The  vapor  portion  of  the  P-F  equation-of-state  is  used  for 
densities  less  than  the  critical  density  when  the  point  lies 
outside  the  liquid-vapor  mixed-phase  region. 


3.2.2  Vapor-Liquid  Mixed-Phase .  Pressure  as  a  function 
of  temperature  in  the  mixed-phase  region  is  obtained  by  using 
the  Clapeyron-Clausius  equation: 


dP  _  AH 

dT  -  T  -  V*) 


(10) 


where  AH  is  the  heat  of  vaporization  at  temperature  T.  The 
subscripts  v  and  Z  refer  to  the  vapor  and  liquid  phases, 
respectively,  at  temperature  T.  The  specific  volume  of  the 
liquid  phase  is  assumed  to  be  negligible  in  comparison  with  that 
of  the  vapor  phase 


Vv  ~  V* 


RT 

P 


and  the  heat  of  vaporization  is  approximated  by 


(11) 
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AH  *=  AHfi  +  AC  (t  -  Tfi) 


(12) 


where  AC  is  the  difference  between  the  specific  heats  of  the 
liquid  and  vapor.  The  subscript  B  refers  to  conditions  at  the 
normal  boiling  point.  With  these  assumptions  the  Clapeyron- 
Clausius  equation  is  intergrated  to  give  pressure  as  a  function 
of  temperature  in  the  liquid-vapor  phase  region. 


H  n 


AHb  -  AC  Tb 
R 


1 

T 


+ 


(13) 


This  expression  is  used  in  the  liquid-vapor  region  for  tempera¬ 
tures  less  than  the  critical  temperature. 


The  critical  pressure  is  not  determined  from  Equation  13  by 
setting  T  =  T  ,  but  is  obtained  by  an  alternate  approach  which 
first  determines  the  critical  compressibility  ratio.  The  critical 
compressibility  ratio  is  determined  by  using  the  form  of  the 
reduced  vapor  pressure  suggested  by  Riedel  and  used  by  Hirschfelder, 
et  al. 


H  n  p  = 


°Rln 


t-0.0838 


(cr-3.75) 


-  35-t  +42  Jtn 


(14) 


where  all  thermodynamic  quantities  are  reduced  variables.  The 
parameter  aR  is  specified  by  requiring  that  Equation  13  satisfy 
the  Clapeyron-Clausius  equation  at  the  normal  boiling  point 


'B 


AH 

RT 


B  1 


(15) 
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The  critical  compressibility  ratio  is  then  obtained  from  the  rela¬ 
tionship  suggested  by  Riedel. 


(16) 


Both  Equation  14  and  Equation  16  which  are  used  to  define  the 
critical  compressibility  ratio  are  based  on  relationships  which 
were  based  on  hydrocarbon  data.  Their  extension  to  metals  does 
seem  questionable;  however,  the  results  of  using  Equation  16  give 
reasonable  agreement  with  the  results  of  more  rigorous  calcula- 
tional  procedures. 

The  line  that  separates  the  mixed-phase  vapor- liquid  region 
from  the  vapor  region  is  constructed  by  setting  Equations  4 
and  14  equal  and  solving  for  the  reduced  vapor  density,  pv,  as 
a  function  of  the  reduced  temperature,  t,  for  t  <  1, 

Pvt/Zc  (l  ”  pv  +  b'  pv2)  -  a*  (t)pv2  -  a(t)pv3 

=  exp  (a-a^)  (1  -  1/t)  +  a ^  l n  t  (17) 

where  =  Ac/R.  Since  pv  is  known  as  a  function  of  t  along  this 

line  and  P  is  given  by  Equation  13,  the  line  that  separates  the 
mixed-phase  vapor-liquid  region  from  the  vapor  region  is  known. 
The  internal  energy  along  this  line  is  found  from  Equation  8. 

For  the  line  chat  separates  the  liquid  from  the  liquid-vapor 
region,  an  expression  proposed  by  Hirschf elder  is  used.  The 
reduced  liquid  density  is 

Pz  it)  =  1  +  c  (1  -  t)1/3  +  d  (1  -  t)  (18) 
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where  c  and  d  are  constants  that  depend  on  the  material.  The 
value  for  the  constants  c  and  d,  which  are  suggested  in  Reference 
5,  are  c  =  1.8  and  d  =  1.75.  These  values  are  appropriate  for 
the  alkali  metals  and  have  been  extended  to  other  metals.  The 
internal  energy  on  the  liquid  side  of  the  vapor-lined  region  is 
determined  from  the  value  of  the  internal  energy  on  the  vapor 
side  (Equation  8)  and  the  Clapeyron-Clausius  equation  (Equation 
10).  Thus,  the  reduced  thermodynamic  quantities  p,  t,  p  and  E, 
are  known  along  the  boundary  of  the  liquid-vapor  region. 

The  critical  density,  pc,  is  determined  by  requiring  that 
the  vapor- liquid  region  join  the  solid  mixed-phase  region  at 
zero  pressure, 

Pc  =  |l  +  .(l-^»*.(x-tJl 

The  critical  pressure  is  then 

Pc  =■  Zc  RTC/VC  (20) 

The  critical  temperature  is  an  input  parameter  of  the  model.  The 
vapor  fraction  F  is  introduced  and  the  equations  are  solved 
iteratively. 


3.2.3  Solid  Phase.  The  equation-of-state  for  the  solid 
phase  in  the  model  is  the  standard  PUFF  solid  equation-of-state: 

P(p,  E)  =  Cp  +  Dp2  +  Sp3  +  (p  +  1)  T  E  p  =  £ - 1  (21) 

po 

The  model  does,  however,  require  that  (3P/3p)£  remain  greater  than 
zero  for  problem  free  application.  This  precludes  using  values  of 
S  which  are  negative.  There  is  also  a  consistency  condition  which 
must  be  satisfied;  it  requires  that  the  pressure  at  the  normal 


boiling  point  be  zero.  The  code  handles  this  internally  by  the 
following  logic: 


a.  If  the  density  at  the  normal  melting  point,  pso,  is 
specified,  the  cede  sets  the  internal  energy  at  the  normal 
boiling  point  (Ego,  so  that  P(pso,  Eso)  =  0. 

b.  If  Eso  is  specified,  then  pso  is  determined  such  that 
P(Pso>  Esq)  =  0* 


3.2.4  Solid-Liquid  Mixed- Phase .  For  the  solid-liquid 
mixed-phase  region  the  Clapeyron-Clausius  equation  is  used 
once  again: 


(22) 


where  AH  is  the  heat  of  fusion.  The  subscripts  l  and  s  refer 
to  the  liquid  and  solid  phase,  respectively.  The  slope  dP/dT 
along  the  melting  line  is  assumed  constant  giving  the  pressure 
as  a  function  of  temperature: 


where  the  subscript  m  refers  to  conditions  at  the  normal  melting 
point.  The  quantities  V^o  and  VSQ  are  the  liquid  and  solid 
density  at  the  zero-pressure  melting  temperature.  As  a  consequence 
of  assuming  that  dP/dT  is  a  constant  along  the  melting  line,  the 
ratio  of  the  energy  difference  to  the  volume  change  between  the 
solid  and  liquid  is  constant. 


AE  . 3(T).-  VT)  *Eo 

AV  =  (T)  -  Vg  (T)  AVq 


(24) 
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The  initial  slope  of  the  internal  energy  as  a  function  of 
temperature  is  known  for  both  the  liquid  and  vapor  phases  since 
the  respective  specific  heats  are  known.  The  average  internal 
energy  of  the  liquid  and  solid  is  assumed  linear  with  temperature 
with  the  proportionality  constant  being  the  average  specific 
heat  of  liquid  and  solid  at  normal  pressure  conditions .  This 
relation  is  given  by 


E(t)  = 


\  (T)  +  Es(T) 
2 


W1W  +  c(t-tJ  (25a) 


where 


ca<V  +  C*'V 


(25b) 


The  form  of  the  expression  for  the  internal  energy  of  the  liquid 
and  solid  phase  is  constructed  so  that  d^/dT  =  and  dEg/dT  =  Cg 
at  temperatures  close  to  T^.  At  large  temperatures  dE^  /dT  *=  C 
and  dE  /dT  88  C.  The  suggested  forms  that  satisfy  these  conditions 
are 


(26c) 
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The  melting  entropy  is  given  by 


AH  ^(T)-Es(T)+P 


vt  (T)  -  Vs  (T) 


for  the  melting  entropy 

] 

AH 

T 

which  reduces  to 


AH 

T 


T 

Equation  23 

provide 

*(E*  •  Es) 

T_  _  x 

T  A 
m 

T 

E.  -  E 
£  s 

T 

m 

(27) 


(23) 


(29) 


Since  AC  is  generally  a  small  quantity,  the  difference  in  internal 
energy  is  nearly  constant.  Therefore,  the  melting  entropy  is 
also  nearly  constant  and  approximately  equal  to 


p  —  P  p  —  P 

AH  _  “£ _ 2s  _  £  o  so 

ip  ij»  ~  ip  W'J  / 

m  m 

The  density  on  the  solid  melting  curve  is  obtained  by 
requiring  that  the  pressure  be  continuous  across  the  solid  melt¬ 
ing  curve.  This  means  that  Equations  21,  and  23  must  be  solved 
simultaneously.  The  density  on  the  liquid  melting  curve  may  then 
be  obtained  by  using  the  Claperyon-Clausius  equation. 

In  the  solid-liquid  phase  the  volume  V  and  internal  energy 
E  are  related  through  the  liquid  fraction  F: 


V  «*  F  Vt 

+  (1  -  F)  Vg 

(31a) 

E  =  FE£ 

+  (1  -  F)  E 

b 

(31b) 
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The  equations  are  solved  iteratively  to  define  P(V,E)  in 
the  solid-liquid  region. 

3.2.5  Liquid  Phase.  The  equation-of-state  formalism  up 
to  this  point  has  determined  all  the  thermodynamic  variables  on 
the  boundary  of  the  liquid  region.  The  equation  of  state  in  the 
liquid  region  is  then  developed  by  using  empirical  curves  along 
which  the  internal  energy  is  constant. 

Two  empirical  expressions  are  suggested.  The  first  assumes 
that  the  pressure  is  a  linear  function  of  density  for  a  given 
energy 


P^P,  E)  = 


P  -  P 
a  m  ZB 

p£m  “  P£B 


PHB  Pfcm  ~  pfcm  PfcB 
p£m  ~  PZB 


where  the  subscripts  m  and  B  refer  to  the  melting  and  boiling 
sides  of  the  liquid  region.  For  E  >  Ec  the  density  becomes 
Pc  and  P^  B  becomes  Pv  (E ,  pc)  where  the  V  refers  to  the  vapor 
regie  ».  The  second  expression  assumes  that  the  logarithm  of  the 
pressure  is  a  linear  function  of  the  logarithm  of  the  density, 

«-  n  f  WP£B) 

i”  p2  (o,  ei  =  irr7rf‘tp  (33> 

^  X,n  (pJlm/pJlB) 

_  *n  pHB*n  PZm  '  *n  p&m&n  P*B 
2,11  (pim/piB) 

The  actual  relation  used  for  the  equation  of  state  in  the 
liquid  region  is  obtained  by  using  a  linear  combination  of  the 
logarithms  of  the  two  pressures,  P^  and 


I 
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an  P  { p ,  E) 


F(EHn  P2  (p,  E)  +  (1  -  F (E) )  £n  P2  (p,  E)  (34) 


where  0  <  F(E)  <1  and  is  determined  by  the  condition 


=WE> 

P|lra(E)  '  plo 


(35) 


By  choosing  F (E)  in  this  manner  P^  is  predominant  at  low  energies, 
while  P2  is  predominant  at  high  energies. 


3.2.6  Discussion  of  the  Model.  The  PHILCO-FORD  equation 
of  state  model  was  developed  for  materials,  such  as  metals, 
which  melt  rather  than  sublime  and  whose  vapor  phase  is  composed 
of  atomic  species.  Since  the  pressure  of  the  triple  point  for 
metals  is  quite  low,  orders  of  magnitude  below  one  atmosphere, 
the  model  has  neglected  the  sublimation  region  of  the  phase 
diagram.  The  thermodynamic  quantities  which  are  measured  at  one 
atmosphere  are  assumed  to  occur  at  zero  pressure. 

The  vapor  equation— of-state  is  based  on  a  modified  van  der 
Waals  equation.  This  model  in  the  vapor  region  should  be  ade¬ 
quate  since  both  the  experimental  data  and  advanced  theoretical 
calculations  have  demonstrated  that  it  is  applicable  in  this 
region.  Electronic  contributions  to  the  specific  heat  have  been 
neglected,  but  at  the  dose  levels  of  interest  for  this  report 
this  contribution  is  negligible.  Also,  the  formalism  does  exist 
for  including  these  considerations  within  the  model. 
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The  envelope  of  the  vapor- liquid  region  is  determined  by 
us.ng  relations  based  on  hydrocarbon  data  and  an  integration  of 
the  Clapeyron-Clausius  equation  with  assumptions  which  are  not 
justified  over  the  total  region  of  integration.  Both  of  these 
objections  make  the  determination  of  the  envelope  questionable. 
However,  the  results  of  the  model  give  values  for  the  critical 
parameters,  which  are  in  reasonable  agreement  with  values 
determined  from  experiments  and  from  advanced  theoretical 
considerations.  This  tends  to  support  the  validity  of  the 
assumptions  used  to  construct  the  envelope  of  the  liquid- vapor 
region  and  the  equation-of-state  model  in  the  liquid- vapor  region. 

The  equation-of-state  in  the  solid  region  is  equivalent 
to  the  PUFF  treatment  for  a  solid,  which  has  withstood  the  test 
of  time  quite  well. 

In  the  mixed-phase  liquid-solid  region  the  melting  entropy 
is  nearly  constant.  This  assumption  has  been  shown  to  be  valid 
for  most  metals  both  by  experiment  and  theoretical  considerations 
(Reference  9) .  Unfortunately  aluminum,  the  metal  of  interest 
in  this  program,  does  not  fall  into  this  class  of  metals 
(Reference  9).  The  specific  heat  of  the  liquid  along  the  melt 
line  is  assumed  to  be  approximately  a  constant. 

The  liquid  equation-of-state  is  the  most  questionable.  The 
relationships  are  entirely  empirical  and  lack  either  experimental 
or  theoretical  data  for  their  justification.  A  much  better 
method  of  treating  both  the  liquid  phase  and  liquid-solid  phase 
would  be  to  incorporate  the  liquid-metal  equation-of-state  based 
on  scaling  proposed  in  Reference  9  and  used  by  Royce  in  the  GRAY 
model,  which  is  discussed  later. 
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In  summary,  the  PHILCO-FORD  equation-of-state  model  has 
provisions  for  treating  the  liquid-solid  and  liquid-vapor  mixed- 
phase  regions  as  well  as  the  standard  solid,  liquid,  and  vapor 
phases.  Its  most  questionable  treatment  is  in  the  liquid-solid 
and  liquid  regions.  The  framework  does  exist  for  modifying 
these  two  regions  with  an  updated  treatment  such  as  suggested 
by  Grover  in  Reference  9. 

3 . 3  GRAY  EOS 

The  GRAY  equation-of-state  model  was  developed  at  the 
Lawrence  Livermore  Laboratory,  Livermore,  California  by  Royce, 
Grover,  Young,  and  Alder  (References  10  through  12) .  The  final 
documentation  is  given  in  Reference  11. 

In  the  solid,  liquid-solid  and  liquid  regimes  the  GRAY 
equation-of-state  model  is  based  on  the  liquid-metal  equation 
of  state  developed  by  Grover  (Reference  9) .  The  vapor  equation 
of  state  has  been  obtained  by  analytically  joining  Grover's 
model  to  the  modified  van  der  Waals  equation-of-state  developed 
by  Young  and  Alder  (Reference  12).  What  follows  is  a  brief 
summary  of  the  GRAY  model. 

3.3.1  Solid,  Solid-Liquid,  and  Liquid  Region.  The  funda¬ 
mental  assumptions  used  by  Grover  in  obtaining  the  equation-of- 
state  in  these  regions  are: 

a.  The  entropy  of  melting  at  constant  pressure  is  a 

constant  independent  temperature  and  pressure. 

b.  The  specific  heat  in  the  liquid  has  the  functional 

form  Cv  =  3R  -  a  T/Tm,  a  =  0.15R,  T  >  T^. 
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c.  The  melting  temperature  as  a  function  of  specific 
volume  can  be  obtained  by  integrating  Lindemann's  relation: 


d  In  T 
_ m 

d  in  v 


2  rs(v> 


were  r  (v)  is  the  volume-dependent  Gruneisen  coefficient 
for  thi  solid  phase . 


In  Reference  9  Grover  carefully  considers  each  assumption  and 
justifies  each  on  the  basis  of  experimental  data,  theoretical 
considerations,  and  the  results  of  numerical  statistical 
mechanics  calculations . 


The  equation  of  state  is  constructed  by  using  assumptions 
(a)  and  (b)  to  obtain  the  Helmholz  free  energy  a.  The  pressure 
and  energy  relations  then  follow  from  the  relations  P  =-(3A/8V)T 
and  E  -  A  +  TS.  The  equation-of-state  is 


A (T,  V)  = 

E(T,  V)  = 

V  P  (T ,  V) 

A(T,  V)  = 
E(T,  V)  = 

V  P(T,  V) 
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where  AT  is  the  temperature  difference  between  the  solidus  tem¬ 
perature  Tg(V)  and  liquidus  temperature  T^(T),Tm  =  (T  +  T^)/2,v 

(T  -  T  (V))/AT(V),  AS  =  entropy  of  melting,  y  =  -  d  in  T  /d  Hn  V, 
a  m 

and  a  =  temperature  coefficient  of  specific  heat.  The  Mie-Grunei 

sen  EOS  for  the  solid  extrapolated  to  temperatures  above  melting 

gives  A  ,  E  ,  and  P  .  The  expression  for  AT  is  obtained  by 
s  s  s 

using  the  Clapeyron-Clausius  equation.  The  integration  of 
Lindemann's  relation  gives  T^,  the  melt  temperature. 

The  specific  heat  decreases  with  increasing  temperature 
until  it  reaches  the  ideal  gas  value  of  3R/2  at  a  temperature 
of  Tg  where  it  remains  constant  for  increasing  temperature.  For 
temperatures  above  T^  the  EOS  takes  on  the  simplified  form 

E(T,  V)  =  Eg(V)  +  |  RT 

V  P(T,  V)  =  V  P_(V)  +  RT 

where 

Vv>  “  Eo(v)  +  Tm(v)  4S  +  l((lf)2  -  X) 

2 

V  PG(V)  =  V  PC(V)  +  MV)  Tm(V)  AS  +  |((|f)  -  1 

and  Eq(V)  is  the  energy  along  the  cold  compression  curve. 


The  GRAY  equation  of  state  has  a  volume-  and  energy-dependent 
Gruneisen  coefficient  y  5  V(8P/8E)^.,  The  solid-phase  Gruneisen 
coefficient  yg  is  a  function  of  volume  only.  Upon  melting  y  ex¬ 
ceeds  y  ;  then,  as  the  temperature  is  increased,  y  slowly  de- 
s 

creases  until  it  reaches  the  ideal  monatomic  gas  value  of  2/3  at 
T  =  Tg  where  T  is  approximately  10  Tm. 

In  summary,  the  GRAY  equation-of-state  model  in  the  solid, 
liquid-solid,  and  solid  region  is  based  on  Grover's  liquid-metal 
EOS  which  is  derived  from  scaling  laws.  The  foundation  of  the 
model  in  these  regions  are  the  scaling  laws  (assumptions  a  through 
c)  which  are  based  on  both  experimental  data  and  theoretical 
considerations . 

3* 3*2  Liquid-Vapor  Region  (Expanded  Region) .  The  liquid- 
vapor  region  uses  a  modified  van  der  Waals  equation-of-state 
developed  by  Young  and  Alder  (Reference  12).  The  Helmholz 
free  energy  is  expanded  about  the  hard-sphere  free  energy  as  a 
perturbation  series.  The  cohesive  energy  is  approximated  by  the 
mean-field  approximation  used  in  the  classical  van  der  Waals 
model.  What  is  obtained  is  a  modified  van  der  Waals  model. 

The  model  has  been  used  to  predict  critical  parameters  of 
the  metallic  elements  and  compared  with  estimates  from  experi¬ 
mental  data  (Reference  12).  The  comparison  with  the  experimental 
data  is  generally  good  in  light  of  the  errors  which  are  present 
in  the  estimates  from  the  experimental  data. 

3.3.3  The  Complete  EOS.  The  complete  EOS  is  obtained  b’ 
analytically  joining  the  Grover  scaling  law  model  to  the  modified 
van  der  Waals  model  of  Young  and  Alder  at  a  specific  volume  of 
1.3  to  1.5  times  the  normal  volume.  The  joining  is  done  in  a 
thermodynamically  consistent  manner.  The  details  and  rationale 
for  the  joining  are  given  in  Reference  11. 
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3.3.4  Discussion  of  the  Model.  The  GRAY  EOS  model  treats 
the  solid,  liquid-solid,  liquid,  and  vapor  regions  of  the  phase 
plane.  The  model  is  based  on  the  results  of  two  models  which 
are  formulated  from  theoretical  considerations  and  which  agree 
reasonably  well  with  existing  experimental  data.  The  model 
does  ignore  the  liquid-vapor  region. 

3.4  COMPARISON  OF  THE  PHILCO-FORD  AND  GRAY  EOS 

Both  the  PHILCO-FORD  and  GRAY  EOS  models  treat  the  solid, 
solid-liquid,  liquid, and  vapor  phases.  The  P-F  model  also 
allows  the  existence  of  the  liquid-vapor  portion  of  the  phase 
plane.  This  section  describes  some  of  the  similarities  and 
differences  of  the  two  models. 

3.4.1  Solid  Region.  The  models  contain  similar  treatments 
of  the  solid  phase.  Both  models  follow  closely  the  Mie-Gruneisen 
equation-of-state  in  the  solid  phase.  The  P-F  model  assumes  that 
the  Gruneisen  coefficient  y  is  constant  while  the  GRAY  model  has 
a  linear  dependence  of  y  on  the  specific  volume. 


3.4.2  Liquid- Solid  Region.  The  models  differ  in  the 
manner  in  which  the  melt  line  is  established.  The  P-F  model 
assumes  that  dP/dT  along  the  melt  line  is  constant  and  equal 


to  the  value  obtained  from  the  Clapeyron-Clausius  equation 
evaluated  at  one  atmosphere.  The  GRAY  model  first  defines  a 


parameter  called  the  melting  temperature 


T  (V) 
m 


Ts  (V)  +  T^  (V) 
2 


where  the  subscripts  s  and  l  refer  to  the  solidus  and  liquidus 
respectively.  The  value  of  T^(V)  is  obtained  by  using  Lindemann's 
law.  An  estimate  of  AT  =  T£(V)  -  Ts (V)  is  obtained  from  the 
Clapeyron-Clausius  equation. 
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The  GRAY  EOS  assumes  that  the  entropy  change  on  melting  is 
a  constant.  The  entropy  of  melting  in  the  P-F  model  is  nearly 
constant  because  of  the  assumptions  made  in  deriving  the  model 
in  the  melt  region.  Both  models  show  an  increase  in  the  Gruneisen 
coefficient  during  melt  for  aluminum. 

3.4.3  Liquid  Region.  The  GRAY  EOS  is  divided  into  two 
subregions.  The  first  is  referred  to  as  the  liquid  region  and 
it  is  for  temperatures  T  <  Tg  (~  10  T^) .  In  this  region  the 
specific  heat  and  Gruneisen  coefficient  are  decreasing  as  the 
temperature  increases.  For  temperatures  above  TQ  the  material 
is  called  a  hot  liquid  with  Cv  =  3R/2  and  y  =  2/3,  the  monatomic 
perfect  gas  value. 

The  P-F  model  uses  empirical  curves  to  span  this  region. 

These  curves  have  no  theoretical  justification.  The  Gruneisen 
coefficient  at  normal  value  remains  above  twr  even  for  doses  of 
as  much  as  20,000  cal/gm.  The  Gruneisen  does,  however,  become 
smaller  for  an  expanded  volume. 

3.4.4  Liquid-Vapor  Region.  The  GRAY  model  does  not 
consider  the  liquid-vapor  region  of  the  phase  plane  as  a  separate 
state.  The  GRAY  model  uses  an  EOS  developed  by  Young  and  Alder 
which  predicts  values  of  the  critical  constants  close  to  the 
experimentally  observed  values .  This  would  seem  to  indicate  that 
GRAY  provides  an  adequate  representation  of  the  EOS  in  the 
critical  region.  However,  the  results  of  the  comparison  of  GRAY, 
PUFF-V  and  P-F  for  aluminum  described  in  the  next  section  do  not 
support  this  conclusion. 
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The  construction  of  the  envelope  of  the  liquid- vapor  region 
in  the  P-F  model  is  not  exact.  However,  the  model  predicts  values 
of  the  critical  constants  which  are  close  to  the  experimentally 
observed  values.  The  model  does  provide  reasonable  definition 
of  the  liquid-vapor  region. 

3.4.5  Vapor  Region.  Both  models  use  a  modified  van  der 
Waal  equation-of-state  in  the  vapor  region.  The  coefficient  of 
attractive  potential,  a,  in  the  van  der  Waal  EOS,  which  was  used 
in  the  numerical  example  for  aluminum  given  m  the  next 
section,  are  comparable  for  each  model. 


aGRAY 


0.06457 


mbar  cm 
2 

gm 


an  =  0.07548 

P-F 


mbar  cm 
gm^ 


However,  the  pressure,  sound  speed,  and  Gruneisen  parameters  do 
not  agree  in  the  expanded  region,  Figures  21  through 2 3  . 


3.4.6  High-Temperature  Behavior.  At  high  temperatures 
(on  the  order  of  2  ev  or  above)  the  electronic  contributions 
to  the  EOS  become  important.  Neither  the  GRAY  nor  the  P-F 
model  takes  into  account  the  electronic  contributions  to  the 
EOS.  Thus  application  of  both  models  to  extremely  high- temperature 
experimental  data  is  questionable. 
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SECTION  4 

COMPARISON  OF  THE  GRAY,  PHILCO-FORD,  AND  PUFF-V 
EQUATION-OF-STATE  MODELS  FOR  ALUMINUM 

4 . 1  INTRODUCTION 

A  comparison  of  the  GRAY,  PHILCO-FORD,  and  PUFF-V  EOS 
models  was  carried  out  for  aluminum.  The  input  parameters  used 
for  each  of  the  EOS  models  are  given  in  Tables  2  through  4. 

4.2  REFERENCE  DENSITY 

The  pressure  predicted  by  each  EOS  model  as  a  function  of 
internal  energy  density  for  aluminum  at  reference  density  is 

shown  in  Figure  14.  The  pressures  predicted  by  the  P-F  model 

\ 

for  aluminum  at  reference  density  are  higher  than  those  pre¬ 
dicted  by  either  the  GRAY  or  PUFF-V  EOS  for  internal  energy 
densities  above  400  cal/gm.  For  internal  energy  densities 
above  3500  cal/gm,  the  P-F  pressures  are  at  least  a  factor  of 
two  above  those  of  GRAY.  The  PUFF-V  and  P-F  models  are  in 
reasonable  agreement  for  energy  densities  less  than  3000  cal/gm. 

In  the  PUFF-V  EOS  model,  aluminum  at  reference  density 
remains  as  solid  as  the  internal  energy  density  is  increased. 

Hie  pressure  as  predicted  by  PUFF-V  is  proportional  to  the 
internal  energy  density  with  the  product  of  the  Gruneisen 
coefficient  times  the  density  being  the  proportionality 
constant  (P  =  p  T  E) . 
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TABLE  2 


PHILCO-FORD  EOS  INPUT  PARAMETERS  FOR  ALUMINUM 


Description 

Units 

Value 

Atomic  weight 

gm/mole 

26.98 

Normal  density 

gm/cm3 

2.71 

Critical  temperature 

K° 

8000.00 

c  coefficient  in  liquid 
density  equation 

None 

1.8 

d  coefficient  in  liquid 
density  equation 

None 

1.75 

kQ  parameter  in  vapor  EOS 

None 

0.0 

Melting  temperature  T^ 

K° 

932.00 

Solid  enthalpy  at  TM 

kca 1/mole 

4.291 

Liquid  enthalpy  at  TM 

kcal/mole 

6.838 

Vapor  enthalpy  at  TM 

kcal/mole 

81.177 

Solid  density  at  TM 

gm/cm3 

2.537 

Liquid  density  at  TM 

gm/cm3 

2.380 

Boiling  temperature  Tfi 

K° 

2736.00 

Liquid  enthalpy  at  Tfi 

kcal/mole 

19.466 

Vapor  enthalpy  at  Tn 

kcal/mole 

90.149 

C  coefficient  in  solid  EOS 

kbar 

720.00 

D  coefficient  in  solid  EOS 

kbar 

1313.00 

S  coefficient  in  solid  EOS 

kbar 

2008.00 

Gruneisen  coefficient  in 

None 

2.18 

solid  EOS 
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TABLE  3 
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TABLE  4 


PUFF-V  EOS  INPUT  PARAMETERS  FOR  ALUMINUM 


Description 

Units 

Value 

Density 

gm/cm'* 

2.71 

C  coefficient  in  solid 

EOS 

mbar 

0.720 

D  coefficient  in  solid 

EOS 

mbar 

1.313 

S  coefficient  in  solid 

EOS 

mbar 

2.008 

Gruneisen  coefficient 

None 

2.13 

K  parameter  in  vapor  EOS 

None 

0.1 

Sublimation  energy 

mbar 

0.3232 

PRESSURE,  mbar 


PHiLCO-FORD  EOS  MODEL 


GRAY  EOS  MODEL 
PUFF  V  EOS  MODEL 


Hot  _ 

_ _ (Gray  t0S 


Energy  density,  cal/gm 


Figure  14  Comparison  of  predicted  pressure  as  a  function  of 
energy  density  for  aluminum  at  reference  density. 
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Both  the  GRAY  and  P-F  EOS  models  allow  aluminum  to  undergo 
a  melting  phase  transition  at  reference  density  for  sufficiently 
high  internal  energy  densities.  GRAY  predicts  the  onset  of 
melting  at  200  cal/gm  with  full  melt  occurring  at  approximately 
300  cal/gm#  while  the  corresponding  values  obtained  from  P-F  are 
300  cal/gm  for  start  melt  and  600  cal/gm  for  end  melt. 

Hie  Gruneisen  coefficient  r  =  V(3P/3E)v  as  a  function  of 
internal  energy  density  for  aluminum  at  reference  density  is 
shown  in  Figure  15.  As  pointed  out  previously,  the  Gruneisen 
coefficient  in  PUFF-V  is  independent  of  internal  energy  and 
equal  to  a  constant. 

Both  the  P-F  and  GRAY  EOS  models  show  an  increase  in  the 
Gruneisen  coefficient  during  the  melting  phase  transformation 
at  constant  volume.  The  Gruneisen  coefficient  during  the  melt 
phase  trans formation  at  constant  volume  is  larger  than  either  the 
solid  or  the  liquid  Gruneisen  coefficient  at  energies  adjacent 
to  the  liquid-solid  region. 

For  energies  above  2500  cal/gm  the  Gruneisen  coefficient  in 
the  GRAY  model  is  equal  to  the  ideal  monatomic  gas  value  of  2/3. 
The  rationale  for  the  decrease  in  T  to  this  limiting  value  is 
given  in  Reference  9. 

The  sharp  increase  in  the  P-F  Gruneisen  coefficient,  which 
occurs  at  3000  cal/gm,  is  due  to  the  manner  in  which  the  model 
is  constructed  in  the  liquid  phase  region.  For  the  liquid  region 
the  pressure  is  determined  by  the  thermodynamic  conditions  on 
the  boundaries  of  the  liquid  phase  region  (see  Figure  16).  For 
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Figure  15  Comparison  of  the  Gruneisen  coefficient  as  a  func¬ 
tion  of  energy  density  for  aluminum  at  reference 
density. 
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Figure  16  P-V  phase  diagram  with  liquid  region  boundaries 
defined  by  a  solid  line. 

internal  energy  densities  below  3000  cal/gm,  pressures  in  the 
liquid  phase  region  are  determined  from  thermodynamic  conditions 
on  the  liquid  side  of  the  liquid-solid  region  (line  A,  Figure 
16)  and  conditions  on  the  envelope  of  the  vapor-liquid  region 
(line  B,  Figure  16).  When  the  internal  energy  densities  are 
above  3000  cal/gm,  the  liquid  phase  pressures  are  determined 
from  thermodynamic  conditions  in  the  vapor  phase  when  p  =  pcritical 


(line  C,  Figure  16)  and  conditions  on  the  boundary  of  the 
liquid-solid  region  (line  A,  Figure  16) .  The  construction 
of  the  EOS  surface  in  the  liquid  region  is  one  of  the  drawbacks 
of  the  P-F  model.  There  is  no  sound  physical  reason  to  expect  a 
discontinuous  change  in  the  pres sure -energy  coupling  coefficient 
for  reference  density  aluminum  at  an  internal  energy  density  of 
3000  cal/gm. 

Another  thermodynamic  quantity  which  influences  the  material 
response  is  the  sound  speed  /(3P/3p) The  sound  speed  predicted 
by  each  model  as  a  function  of  energy  density  for  aluminum  at 
reference  density  is  shown  in  Figure  17.  The  comparison  of 
sound  speeds  is  similar  to  the  Gruneisen  coefficient  comparison. 
For  a  given  value  of  energy  density  the  model  that  predicts  the 
highest  sound  speed  also  predicts  the  highest  Gruneisen  co¬ 
efficient,  while  the  model  that  predicts  the  lowest  sound  speed 
also  predicts  the  lowest  coefficient. 

The  P-F  model  predicts  a  decrease  in  sound  speed  as 
aluminum  melts  at  reference  density.  The  sound  speed  then 
increases  as  aluminum  is  heated  through  the  melt  region.  At 
the  end  of  the  melting  phase  transformation  at  reference  density, 
the  sound  speed  in  the  liquid  phase  is  higher  than  in  the  melt 
phase.  Physically  this  means  that  the  slope  of  the  release 
adiabat  —|  =  - (pc)  is  discontinuous  at  both  the  begin  melt 

and  end  meltf^a^* 


In  the  GRAY  EOS  model  the  sound  speed  is  continuous  across 

the  melting  phase  transformation  and  thus  the  slope  of  the 

3P  2 

release  adiabat  -^|  =  -(pc)  is  continuous  at  both  begin 

melt  and  end  melt.  a<^ak* 
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Figure  17  Comparison  of  predicted  sound  speed  as  a  function  of 
energy-density  for  aluminum  at  reference  density. 


It  is  apparent  from  Figures  14  through  17  that  the  GRAY, 
PHILCO-FORD,  and  PUFF-V  ECS  models  are  in  significant  disagree¬ 
ment  in  the  predicted  values  of  pressure,  Gruneisen  coefficient, 
and  sound  speed  as  a  function  of  energy  density  for  reference 
density  aluminum.  The  question  which  arises  is  the  following: 
can  the  data  (total  impulse  and  propagated  stress  pulse  measure¬ 
ments)  be  used  to  resolve  the  differences  and  choose  which  EOS 
model  is  more  appropriate,  or  determine  a  new  set  of  more  appro¬ 
priate  input  parameters? 


In  order  to  answer  the  above  question,  consider  the  case  of 
melt-dominated  material  response.  Assume  an  EOS  of  the  form 

p  -  f  (y  ♦  p r « 

where  P  is  the  pressure,  P  is  the  density,  Pq  the  reference 
density,  r  is  the  Gruneisen-energy  coefficient,  and  E  is  the 
internal  energy  density.  For  the  linearized  problem 


where  i>  is  the  fluence,  t  the  deposition  time,  c  the  sound  speed, 
and  E(s)  is  the  normalized  energy  deposition  profile.  When 
Equation  35  is  specialized  to  a  linear  deposition  profile  with 
range  r 
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where  Eq  =4>Eo  is  the  maximum  deposited  energy.  If  we  make  the 
additional  assumption  that  the  material  is  removed  to  half  range, 
then  the  impulse  has  been  shown  (Reference  2)  to  be  given  by 


I 


_r_ 

2C 


CT 


(37) 


In  obtaining  both  Equation  35  and  37  the  sound  speed  c  and 
Gruneisen  coefficient  r  were  assumed  to  be  a  constant.  In 
applying  these  relations  to  a  case  where  both  r  and  c  are  not 
constant,  a  more  appropriate  parameter  to  use  is  the  effective 
or  average  value  of  sound  speed  and  Gruneisen  coefficient. 

These  relations  (Equation  36  and  Equation  37)  do  provide  a 
useful  qualitative  way  of  interpreting  the  effect  of  the  sound 
speed  and  Gruneisen  coefficient  on  the  peak  propagated  stress 
and  impulse  for  melt-dominated  material  response. 

For  the  data  on  aluminum  taken  as  part  of  this  program,  the 
deposition  time  (t)  ,  range  (r)  ,  and  fluence  (4>)  are  known  to 
within  the  experimental  error.  The  thermodynamic  parameters 
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which  to  first  order  effect  the  material  response  are  the 
Gruneise-  coefficient  r and  sound  speed  c.  From  Equation  36 
and  Equation  37  it  is  seen  that  for  a  given  deposition  time  an 
EOS  model  that  predicts  a  higher  value  of  both  r  and  c  may  give 
approximately  the  same  value  for  peak  propagated  pressure 
and  impulse  as  a  model  that  predicts  a  lower  value  of  both 
I'  and  c . 

For  aluminum  the  PHILCO-FORD  EOS  generally  predicts  the 
highest  value  of  Gruneisen  coefficient  (Figure  15)  and 
sound  speed  (Figure  17) ,  while  the  GRAtf  EOS  model  generally 
p.edicts  the  lowest  value  of  both  r  and  c.  In  general  the  PUFF-V 
EOS  model  gives  intermediate  values  of  r  and  c.  A  larger  value  of  r 
alone  tends  to  give  a  larger  predicted  value  of  peak  propagating 
stress  (Equation  36)  and  impulse  (Equation  37) .  However,  j. 
larger  value  of  c  tends  to  give  a  lower  predicted  value  of 
peak  propagating  stress  and  impulse.  The  two  effects  tend  to 
compensate  since  a  higher  Gruneisen  coefficient  combined  with 
a  higher  sound  speed  (Equation  36  and  Equation  37)  may  give 
approximately  the  same  value  of  peak  propagating  stress  and 
total  impulse  as  a  lower  value  of  both  T  and  c. 

Thus  we  see  that  from  measurements  of  peak  propagated  stress 
and  total  impulse  alone  it  is  difficult  to  resolve  the  differences 
between  the  EOS  models  and  choose  which  EOS  model  is  more  appro¬ 
priate.  It  may  be  pointed  out  that  Equation  36  and  Equation  37 
are  most  applicable  to  melt-dominated  material  response.  The 
differences  in  the  predictions  of  each  model  for  the  high  fluence 
data  are  discussed  in  the  next  section.  As  the  peak  stress  propa¬ 
gates  it  will  be  attenuated,  due  to  hydrodynamic  attenuation. 

This  attenuation  tends  to  bring  the  predictions  of  the  EOS  models 
closer  together.  In  other  words,  the  model  that  predicts  the 
higher  peak  propagating  stress  generally  predicts  a  larger  hydro- 
dynamic  attenuation  near  the  front  surface. 
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4.3  CRITICAL  DENSITY 


A  comparison  of  the  GRAY,  PHILCO-FORD,  and  PUFF-V  EOS  models 

3 

was  also  carried  out  for  aluminum  at  a  density  of  0.6  gm/cm 
which  is  near  the  critical  density.  The  results  of  the  com¬ 
parison  are  shown  in  Figures  18  through  20. 

Both  the  GRAY  and  PUFF-V  EOS  models  only  allow  for  the  vapor 
phase  to  exist  at  this  density.  The  PHILCO-FORD  EOS  model  does 
include  the  vapor-liquid  mixed  phase  region  of  the  phase  plane 
and  allows  the  material  to  be  a  liquid  or  a  mixed  phase  (vapor- 

3 

liquid)  at  this  density  (0.6  gm/cm  ).  In  the  PHILCO-FORD  model 
the  dividing  line  between  a  liquid  and  vapor  is  p  =  Pcritica3 
(see  Figure  16) .  The  critical  density  for  aluminum  has  not 
been  experimentally  determined  but  is  on  the  order  of  0.5  to  0.7 
gm/cnf*.  The  PHILCO-FORD  EOS  model  predicts  a  critical  density 

3 

of  0.56  gm/cm  . 

For  a  given  energy  density  the  PHILCO-FORD  EOS  model 
predicts  a  larger  value  of  pressure  (Figure  18) ,  Gruneisen  co¬ 
efficient  (Figure  19) ,  and  sound  speed  (Figure  20)  than  either 
the  GRAY  or  PUFF-V  EOS  models.  Of  the  three  models  the  GRAY  EOS 
model  always  predicts  the  smallest  value  for  the  three  quantities 
compared.  The  PUFF-V  EOS  model  predicts  values  that  are  inter¬ 
mediate  to  the  other  two  models. 


For  energy  densities  above  3000  cal/gm  the  PUFF-V  and 
PHILCO-FORD  EOS  models  are  almost  equivalent.  Not  only  do  the 
models  predict  approximately  the  same  value  of  pressure  and 

pressure-energy  coupling  coefficient  for  the  same  internal  energy 

(3P 
3v 

are  aiso  approximately  one  same. 


adiab. 


- (pc) 
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Figure  18  Comparison  of  predicted  pressure  as  a  function  of 
energy  density  for  aluminum  at  a  density  of 
o  =  0,6  gram/ cm^. 
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Figure  20  Comparison  of  predicted  sound  speed  as  a  function  of 
energy  density  for  aluminum  at  a  density  of 

p  =  0.6  gm/cnP. 
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The  GRAY  EOS  model  predicts  significantly  lower  pressures 
than  either  the  PUFF-V  or  PHILCO-FORD  EOS  models. 

4.4  EXPANDED  DENSITY 

A  comparison  of  the  three  EOS  models  was  carried  out  for 
aluminum  with  a  density  of  0.3  gm/cm  .  This  density  represents 
approximately  a  9  to  1  expansion  of  aluminum  from  its  reference 
density  condition.  The  results  of  the  comparison  are  summarized 
in  Figures  21  through  23.  The  comparison  follows  the  general 

3 

trend  of  the  two  previous  comparisons  at  a  density  of  0.6  gm/cm 

3 

and  2.71  gm/cm  .  For  a  given  energy  density  the  PHILCO-FORD  EOS 
model  predicts  the  largest  value  of  pressure,  Gruneisen  co¬ 
efficient,  and  sound  speed  while,  of  the  three  models,  GRAY 
once  again  predicts  the  lowest  value  for  the  three  quantities 
investigated. 

For  energy  densities  above  3200  cal/gm  the  predictions  of 
the  PUFF-V  and  PHILCO-FORD  models  are  almost  i^ntical.  For  a 
given  energy  uensity  greater  than  3200  cal/gm,  the  pressure  and 
pressure-energy  coupling  coefficient  predicted  by  each  model  are 
approximately  the  same.  Since  the  models  also  predict  identical 
sound  speeds,  the  slopes  of  the  release  adiabats  |^| 
are  also  the  same  for  both  of  the  models.  V 

The  only  difference  between  the  two  models  is  in  the  mixed  phase 
liquid-vapor  region  of  the  phase  plane. 

4 . 5  SUMMARY 

In  general,  for  a  given  density  and  internal  energy  density, 
the  PHILCO-FORD  EOS  model  predicted  higher  values  of  pressure, 
Gruneisen  coefficient,  and  sound  speed  than  either  the  PUFF-V 
or  CRAY  EOS  models.  The  GRAY  EOS  predicts  lower  values  of  these 
quantities  for  almost  all  of  the  range  of  densities  and  internal 
energy  densities  investigated. 
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Figure  21  Comparison  of  predicted  pressure  as  a  function  of 
energy  density  for  aluminum  at  a  density  of 
p  =  0.3  gm/cnw. 
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Figure  23  Comparison  of  predicted  sound  speed  as  a  function  of 
energy  density  for  aluminum  at  a  density  of 
p  =  0.3  gm/cm3. 


By  examining  Figures  18  through  20  and  Figures  21 
through  23,  it  can  be  concluded  that  in  the  case  of  aluminum 

3 

for  densities  less  than  0.6  gm/cm  and  energies  greater  than 
3200  cal/gm  the  PUFF-V  EOS  model  and  the  PHILCO-FORD  EOS  model 
are  approximately  equivalent. 

The  comparison  of  the  three  EOS  models  did  point  out  some 
significant  disagreement  in  the  predicted  values  of  pressure, 
Gruneisen  coefficient,  and  sound  speed,  especially  for  reference 
density  conditions.  It  was  concluded  that  it  is  difficult  to 
resolve  the  differences  between  the  EOS  models  and  choose  which 
EOS  model  is  more  appropriate  from  measurements  of  peak  pro¬ 
pagated  stress  and  total  impulse  alone.  The  differences  in  the 
models  tend  to  be  self-compensating  (see  Section  4.2).  It  is 
now  possible  to  design  special  experiments  to  resolve  some  of 
the  differences. 
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SECTION  5 


COMPARISON  OF  THE  EXPERIMENTAL 
DATA  WITH  EXISTING  EQUATION-OF-STATE  MODELS 

The  objective  of  this  program  was  to  investigate  the 
radiation-induced  response  of  aluminum  for  cases  where  a  signi¬ 
ficant  portion  of  the  blowoff  material  was  in  the  mixed 
liquid-vapor  phase.  Figure  24  gives  the  deposited  energy  as  a 
function  of  depth  in  aluminum  for  the  highest  fluence  condition 
obtained  as  part  of  this  program.  Overlayed  on  this  deposition 
profile  are  the  depths  of  the  various  phase  regions  as 
determined  from  the  atmospheric  values  of  enthalpy  for  begin 
melt,  complete  melt,  begin  vaporization,  and  complete 
vaporization.  As  can  be  seen  from  Figure  24,  the  amount  of 
material  that  passes  through  the  mixed  liquid-vapor  phase  is  a 
large  fraction  of  the  material  that  is  removed;  (removal  depth 
is  approximately  equal  to  the  melt  depth) .  It  must  be  noted 
that  Figure  24  uses  the  atmospheric  values  of  the  enthalpy  for 
the  phase  transformation,  and  instantaneous  deposition.  The 
actual  experimental  results  are  for  a  finite  deposition  time 
and  the  pressures  in  the  heated  region  are  greater  than  one 
atmosphere  when  impulse  is  being  generated.  However,  if  the 
results  of  the  numerical  calculations  using  the  PHILCO-FORD  EOS 
model  (the  only  model  studied  that  contained  the  liquid-vapor 
region)  are  examined,  we  find  that  Figure  24  provides  a  conser¬ 
vative  estimate  of  the  amount  of  material  which  passes  through 
the  mixed  liquid-vapor  phase  region. 


Vapor 


For  the  conditions  given  in  Figure  24  approximately  60 

percent  of  the  blowoff  material  is  in  the  liquid- vapor  state. 

2 

If  the  fluence  is  dropped  to  100  cal/cm  and  the  mean  electron 
energy  and  average  angle  remains  the  same  as  Figure  24; 
approximately  60  percent  of  the  blowoff  material  is  still  in  the 

liquid-vapor  state.  This  percentage  drops  until  for  a  fluence 

2 

of  42  cal/cm  only  a  small  fraction  of  the  blowoff  material  is 
in  the  liquid-vapor  state. 

Thus  for  the  experimental  data  taken  at  fluences  above 
2 

100  cal/cm  ,  significant  amounts  of  the  blowoff  material  are  in 
the  liquid-vapor  phase  region.  For  the  condition  given  in 
Figure  24  the  peak  pressure  in  the  heated  region  is  on  the  order 
of  260  kbar  and  this  pressure  occurs  at  a  depth  of  approximately 
0.017  cm.  The  maximum  pressure  of  aluminum  in  the  liquid-vapor 
state  is  on  the  order  of  4  kbar.  Thus  w-  see  that  the  material 
which  nad  the  peak  pressure  of  260  kbar  must  expand  until  its 
pressure  has  dropped  to  at  least  1  kbar  before  it  enters  the 
liquid-vapor  state.  Since  4  kbar  is  only  1.5  percent  of  260 
kbar  there  is  some  question  of  the  importance  of  the  mixed 
liquid-vapor  phase  in  determining  the  material  response.  For 

2 

the  experimental  data  taken  at  fluence  on  the  order  of  50  cal/cm 
the  peak  pressure  in  the  heated  region  is  approximately  80  kbar. 
Even  at  this  lower  fluence  the  maximum  pressure  of  the  liquid- 
vapor  state  is  only  5  percent  of  the  pressure  due  to  direct 
energy  deposition.  Thus  we  see  that  the  material  response  of 
aluminum  for  energy  loading  conditions  for  which  a  significant 
portion  of  the  blowoff  is  in  the  liquid-vapor  state  at  relatively 
low  pressures  may  not  be  dominated  by  mixed  phase  (liquid-vapor) 
effects,  since  the  maximum  pressure  of  the  liquid-vapor  state 
is  only  a  few  percent  of  the  maximum  pressure  due  to  direct  energy 
deposition. 
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The  impulse  data  are  presented  in  Figure  25  as  a  function 
of  fluence.  Also  included  in  Figure  25  is  the  linear  fluid 
model  expression  for  impulse.  The  expression  is  included  in 
Figure  25  to  remind  the  reader  that  each  impulse  measurement 
has  its  individual  set  of  loading  conditions  (e.g.,  mean  electron 
energy,  average  angle  of  incidence,  fluence,  and  deposition 
time)  on  which  impulse  generation  may  depend.  Therefore,  a  plot 
of  impulse  versus  fluence  for  the  electron  beam  data  does  not 
give  a  valid  representation  of  the  scatter  in  the  data,  and  such 
graphs  can  only  be  used  to  present  general  trends. 

The  mean  electron  energy  of  the  flagged  data  point  in 
Figure  25  was  620  keV  and  thus  should  not  be  considered  with 
the  other  data  points . 

A  comparison  of  model  predictions  with  impulse  data, 
generated  in  electron  beams,  must  be  done  on  a  point-by-point 
basis  since  each  impulse  measurement  has,  as  mentioned  before, 
its  individual  set  of  loading  conditions.  From  the  impulse  data 
shown  in  Figure  25  several  representative  data  points  were 
selected  for  a  shot-by-shot  comparison  with  the  theoretical 
predictions  using  the  PUFF-V,  LINEAR-FLUID,  GRAY,  and  PHILCO-FORD 
EOS  models.  The  data  points  selected  are  the  filled-in  points 
in  Figure  25  .  The  selection  of  a  particular  data  point  was 

based  on  the  confidence  of  the  experimentalist  in  the  data  point 
and  in  obtaining  data  points  that  spanned  the  range  of  peak 
doses  for  which  data  existed. 

The  machine  diagnostics  were  used,  as  described  in  Refer¬ 
ence  2,  to  obtain  the  energy  spectrum,  mean  angle  of  incidence 
of  the  electrons,  fluence,  and  effective  deposition  time.  A 
Monte  Carlo  electron  transport  code  was  then  used  to  calculate 
the  deposition  profile  which  was  used  in  the  hydrodynamic  code 
calculations  with  various  EOS  models. 
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The  experimental  data  were  compared  with  the  prediction 
which  used  the  PUFF-V,  LINEAR-FLUID,  GRAY,  and  PHILCO-FORD  EOS 
models.  For  the  theoretical  predictions  each  of  the  EOS  models 
was  combined  with  a  one -dimensional  Lagrangian  finite-difference 
code  (PISCES  1DL) . 

The  LINEAR-FLUID  model  is  given  in  Reference  2  and 
summarized  in  Table  5. 

TABLE  5 

LINEAR-FLUID  CONSTITUTIVE  RELATIONS  FOR  ALUMINUM 
P  =  0.75p  +  0. 405y2  +  (2  +  2p)E,  p  >  0 

P  =  0.75p  +  (2  +  2p)E,  p  <  0 

where 


The  input  parameters  for  the  PUFF-V,  GRAY,  and  PHILCO-FORD  EOS 
model  are  given  in  Tables  2,  3,  and  4  respectively.  A 
von  Mises  vield  model  was  used  with  Y^  =  2.78  kbar.  The  shear 
modulus  was  277  kbar  and  the  spall  strength  was  -13  kbar.  Both 
the  shear  modulus  and  spall  strength  decreased  linearly  with 
internal  energy  density  and  approached  zero  as  the  internal 
energy  reached  the  enthalpy  required  to  reach  the  melting  point 
at  one  atmosphere. 
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The  phase  regions  described  by  each  EOS  model  are  given  in 
Table  6.  It  is  noted  that  only  the  PHILCO-FORD  EOS  model 
attempts  to  describe  all  five  phase  regions  (solid,  solid- liquid, 
liquid,  liquid-vapor,  and  vapor).  The  results  of  the  comparison 
of  the  impulse  data  with  the  predictions  of  the  four  EOS  models 
are  summarized  in  Table  7  and  Figure  26. 

In  Figure  26  the  error  bars  for  each  data  point  have  been 

retained.  As  the  fluence  increases  the  uncertainty  in  the  data 

also  increases.  The  main  source  of  uncertainty  in  the  high- 

fluence  data  is  the  uncertainty  in  the  irradiated  area.  For  the 

electron  beam  tests  impulse  and  fluence  are  not  measured  directly. 

What  is  measured  is  the  total  energy  in  the  beam  (H)  and  the 

total  momentum  imparted  to  the  sample  (AMV)  by  the  material  that 

is  removed.  The  fluence  and  impulse  are  then  determined  by 

dividing  by  the  irradiated  area  A  (i.e.,  <J>  =  H/A  and  I  =  AMV/A) . 

For  a  given  mean  electron  energy  the  total  electron-beam-machine 

output  in  calories  (H)  remains  relatively  fixed.  The  fluence  is 

adjusted  by  using  the  area  which  gives  the  desired  fluence, 

<p  (A  =  H/<J>).  The  data,  which  were  taken  at  higher  and  higher 

fluences,  are  obtained  by  continuously  reducing  the  area.  At 

the  high-fluence  levels  a  relatively  small  uncertainty  in  the 

absolute  value  of  the  irradiated  area  can  give  an  extremely 

large  percentage  uncertainty  in  impulse  and  fluence.  The 

dependence  of  impulse  and  fluence  uncertainty  on  the  uncertainty 

in  area  of  irradiation  (uncertainty  in  crater  radius)  is  shown 

in  Figure  27.  The  area  of  irradiation  for  the  high-fluence 
o 

data  was  <  1  cm  .  From  Figure  27  we  see  that  an  error  of  only 
1  mm  in  determining  the  crater  radius  leads  to  an  uncertainty  of 
approximately  35  percent  in  the  fluence  and  impulse.  The  obvious 
solution  to  the  problem  is  to  use  larger  areas  of  irradiation 
while  maintaining  the  same  fluence.  With  the  installation  of 
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TABLE  6 


PHASE  REGIONS  DESCRIBED  BY  EOS  RESPONSE  MODELS 


Solid- 

Solid 

Liquid 

Liquid 

PUFF  X 

PHILCO-FORD  X 

X 

X 

GRAY  X 

X 

X 

LINEAR  FLUID  X 

X 

Liquid- 

.  Vapor  Vapor 

X 

X  X 

X 
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Fluence,  cal/cnr 


Figure  26  Comparison  of  predictions  of  various  EOS  models  with 
experimentally  observed  impulse  in  aluminum  for  peak 
doses  from  300  cal/gram  to  3800  eal/gram. 


an  external-magnetic  field  beam-control  technique  on  the  Pulserad 
Model  738  (the  electron  beam  machine  on  which  the  data  for  this 
report  were  taken) ,  the  desired  larger  areas  have  become  avail¬ 
able.  The  lines  in  Figure  26  are  only  meant  to  point  out  general 
trends  in  the  comparison.  As  can  be  seen  from  Figure  26  the 
PHILCO-FORD  equation-of-state  model  gives  the  best  representation 
of  the  data. 

In  Figure  26  the  predictions  of  the  PUFF-V  and  LINEAR-FLUID 
ECS  have  been  combined,  since  the  impulse  predictions  of  the  two 
models  for  each  data  point  examined  are  almost  identical  (Table 
7).  As  shown  in  Table  6  the  LINEAR-FLUID  model  does  not 
contain  provisions  for  treating  the  vapor  expansion  as  does  the 
PUFF-V  model.  Tne  equivalence  of  the  predictions  of  the  two 
models  for  peak  doses  less  than  2600  cal/gm  is  not  surprising 
since  it  was  shown  in  Reference  2  that  the  impulse  is  expected 
to  be  melt-dominated  for  peak  dose  up  to  2600  cal/gm.  However, 
the  PUFF-V  EOS  model,  which  includes  a  vapor  state,  and  the 
LINEAR-FLUID  EOS  model,  which  does  not,  predict  approximately 
the  same  impulse  for  peak  dose  up  to  at  least  3800  cal/gm.  Thus 
we  see  the  first  evidence  that  even  for  a  peak  deposited  energy 
of  3800  cal/gm,  which  is  above  the  enthalpy  needed  to  completely 
vaporize  aluminum  (3341  cal/gm)  at  one  atmosphere,  the  impulse 
is  not  dominated  by  liquid-vapor  effects.  Instead,  the  impulse 
production  seems  to  be  dominated  by  the  behavior  of  the  EOS  of 
aluminum  in  the  liquid  region  of  the  phase  plane  (Figure  16) . 

For  a  mean  electron  energy  on  the  order  of  300  keV  and 
fluences  of  less  than  100  cal/cmfc  the  GRAY  and  PHILCO-FORD  EOS 
models  predict  approximately  the  same  value  of  impulse  for  the 

2 

same  energy  loading  conditions.  Above  a  fluence  of  100  cal/cm 
the  predictions  of  the  GRAY  model  are  below  those  of  the  PHILCO- 
FORD  model.  The  difference  in  the  predictions  of  the  two  models 
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becomes  larger  as  the  fluence  increases.  For  the  highest  fluence 

2 

condition  considered  (211  cal/cm  )  the  prediction  of  GRAY  was 

19  percent  below  that  of  PHILCO-FORD.  The  difference  in  the 

impulse  predictions  of  the  two  models  tor  fluencos  above  100 
2 

cal/cm  is  probably  due  to  the  large  differences  between  these 
two  EOS  models  in  the  liquid  region. 

To  investigate  the  reasons  for  the  differences  in  predicted 
impulse  between  the  two  models,  release  adiabats  from  reference 
density  were  obtained  for  aluminum.  The  release  adiabats  from 
reference  density  with  initial  energy  densities  of  300  cal/gm, 
1000  cal/gm,  and  2000  cal/gm  are  shown  in  Figures  28,  29,  and 
30.  The  jump-off  velocity  shown  in  each  figure  is  indicative 
of  the  impulse  produced  by  each  model  as  the  material  expands. 

The  jump-off  velocity  is  given  by 

UJO  “  /  (pc)® 

o  adiab. 

From  Figures  28,  29,  and  30  it  is  apparent  why  the  two  models 
disagree  in  the  value  of  predicted  impulse  for  peak  deposited 
energies  greater  than  1600  cal/gm.  The  jump-off  velocity  and 
the  area  under  the  p-v  curve  (energy  converted  to  kinetic  energy 
of  motion)  are  approximately  equal  for  the  300  cal/gm  and  1000 
cal/gm  release  ac^iabat.  However,  for  the  2000  cal/gm  release 
adiabat,  the  jump-off  velocity  and  the  area  under  the  p-v  curve 
for  the  PHILCO-FORD  model  are  significantly  larger  than  the  jump- 
off  velocity  and  the  area  under  the  p-v  curve  for  the  GRAY  model. 
Thus,  at  an  energy  density  somewhere  between  1000  cal/gm  and 
2000  cal/gm  the  predictions  of  the  two  models  should  start 
disagreeing.  The  wumerical  comparison  of  the  models  places  the 
onset  of  this  disagreement  to  between  a  peak  deposited  energy  of 
1600  cal/gm  and  2500  cal/gm  for  a  mean  electron  energy  on  the 
order  of  300  keV. 
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Figure  29  Comparison  of  release  adiabat  of  aluminum  from 

reference  density  and  an  initial  energy  density  of 
1000  cal/gm. 
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Figure  30  Comparison  of  release  adiabat  of  aluminum  from 

reference  density  and  an  initial  energy  density  of 
2000  cal/gm. 


The  PUFF-V  EOS  model  predicts  lower  values  of  impulse  than 

the  PHILCO-FORD  model  for  all  the  conditions  considered.  For 

the  energy  loading  conditions  obtained  in  this  program  the  PUFF-V 

EOS  model  predicts  lower  values  of  impulse  than  the  GRAY  EOS 

2 

model  for  fluences  less  than  155  cal/cm  .  For  fluences  above 
2 

155  cal/cm  GRAY  predicts  lower  values  of  impulse  than  PUFF-V. 

The  agreement  of  the  predicted  impulse  for  the  LINEAR-FLUID 
and  PUFF-V  EOS  models  for  all  the  data  considered  was  the  first 
indication  that  the  impulse  is  not  dominated  by  liquid-vapor  or 
vapor  effects  for  deposited  energies  up  to  3800  cal/gm.  Another 
indication  was  that  for  all  the  cases  considered  numerically, 
impulse  production  ceased  when  the  material  spalled.  This  was 
independent  of  the  EOS  which  was  used  when  the  nominally  recom¬ 
mended  input  values  for  the  EOS  models  were  used  (Table  2 
through  5 ) . 

The  largest  impulse  delivery  time,  according  to  the  numerical 

calculations,  was  117  nanoseconds  (approximately  three  times  the 

deposition  time)  and  was  obtained  with  the  PHILCO-FORD  model  for 
2 

the  211  cal/cm  (peak  dose  3800  cal/gm)  data  condition.  The 
impulse  delivery  time  for  materials  such  as  SC- 100 8  which  have 
vapor-dominated  impulse  for  comparable  energy  loading  conditions 
is  on  the  order  of  a  few  microseconds  (Reference  14). 

In  order  to  investigate  which  region  of  the  phase  dominates 

2 

the  material  response  behavior,  the  211  cal/cm  data  condition 
was  examined  in  more  detail  with  the  PHILCO-FORD  EOS  model.  The 
peak  deposited  energy  for  this  condition  is  approximately  3800 
cal/gm.  Figure  31  shows  the  release  adiabat  from  reference 
density  for  aluminum  with  an  initial  energy  density  of  4000  cal/gm 
using  the  PHILCO-FORD  model.  As  the  material  expands  it  does 
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Figure  31  Release  adiabat  of  aluminum  from  reference  density 
and  an  initial  energy  density  of  4000  cal/gm  using 
the  Philco-Ford  EOS  model. 


enter  the  liquid-vapor  state  but  at  an  energy  density  of  approxi¬ 
mately  2000  cal/gm  and  a  pressure  of  1.3  kbar.  Thus  we  see  that 
the  liquid-vapor  region  of  the  phase  plane  does  not  dominate  the 
release  adiabat  behavior.  Instead  the  release  adiabat  behavior 
is  dominated  by  the  behavior  of  the  material  in  the  liquid  phase . 
As  used  in  this  report  the  term  liquid  phase  denotes  the  region 
of  the  phase  plane  denoted  as  liquid  in  Figure  16.  At  the 
extreme  pressures  and  energy  densities  in  which  we  are  interested, 
the  terms  liquid  and  vapor  lose  their  significance.  It  is  more 
appropriate  to  refer  to  a  region  of  the  phase  plane. 

2 

The  hydrodynamic  code  calculation  for  the  211  cal/cm 
condition  using  the  PHILCO-FORD  model  provides  additional 
information  on  the  importance  of  the  liquid-vapor  state.  Accord¬ 
ing  to  the  calculation  at  a  time  of  51  nanoseconds  after  irradia¬ 
tion  begins  (deposition  time  is  46  nanoseconds)  all  the  material 
except  for  the  first  zone  in  the  calculation  is  in  the  liquid, 
liquid-solid,  or  solid  state.  The  delivered  impulse  is 
7.63  kilotaps.  This  is  to  be  compared  with  the  total  impulse  of 
9.29  kilotaps  for  an  infinite  time.  Thus  we  see  that  at  least 
82  percent  of  the  impulse  is  delivered  while  the  material  is  in 
the  liquid,  liquid-solid,  or  solid  state.  The  maximum  internal 
energy  density  from  the  calculation  is  approximately  3100  cal/gm 
since  the  material  near  the  front  surface  expands  significantly 
during  the  deposition  time.  Thus  we  see  that,  according  to  the 
PHILCO-FORD  model,  the  impulse  production  seems  to  be  dominated 
by  the  behavior  of  the  EOS  of  aluminum  in  the  liquid  region  of 
the  phase  plane  (Figure  16). 
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The  GRAY  EOS  model  contains  a  parameter  a^,  the  coefficient 
of  attractive  potential  for  the  vapor,  which  can  be  adjusted  to 
evaluate  the  importance  of  the  vapor  and  of  having  a  good 
equation  of  state  in  critical  region.*  As  suggested  by  Royce, 

a  in  GRAY  was  reduced  by  a  factor  of  3.  The  reduced  a  was 

•>  2  * 
used  to  calculate  the  response  of  the  211  cal/cm  condition. 

The  total  impulse  increased  from  7.59  kilotaps  to  11.6  kilotaps 
which  represents  an  increase  of  53  percent  in  impulse  by  decreas¬ 
ing  ay  by  a  factor  of  three.  This  seems  to  indicate  that  a  good 
equation  is  required  in  the  critical  region.  However,  it  must 
be  kept  in  mind  that  the  definition  of  a  vapor  in  GRAY  is  the 

condition  that  the  specific  volume  V  >  1.32  V  .  Also  GRAY  with 

o 

the  nominal  value  of  predicts  much  lower  pressures  than  the 

other  models  (Figures  14,  18,  and  21).  Another  interesting 

2 

phenomenon  occurs  for  the  211  cal/cm  condition  when  ay  is 
reduced  by  a  factor  of  three  in  the  GRAY  EOS.  A  double  wave 
structure  develops  in  the  propagated  stress  pulse.  The  develop¬ 
ment  of  the  double  wave  structure  is  shown  in  Figure  32.  The 
upper  portion  of  Figure  32  is  for  a  time  of  51  nanoseconds  after 
the  start  of  irradiation.  The  normal  compressive  pulse  is 
followed  by  a  strong  rarefaction  which  is  followed  by  another 
compression-rarefaction  wave.  The  first  strong  rarefaction 
moves  into  the  body  and  causes  a  spall  to  develop  in  the 
liquid-solid  phase;  (note  that  the  spall  strength  was  set  equal 
to  zero  for  material  in  the  liquid-solid  phase) .  The  second 
compressive  pulse  recompresses  the  spalled  material  and  then 
propagates  into  the  sample  (lower  portion  of  Figure  32). 

The  question  which  immediately  arises  is  the  following:  Is 
there  any  experimental  evidence  of  the  existence  of  the  double 
wave  structure  discussed  above?  Several  laser  interferometer 
measurements  were  made  with  thin  aluminum  samples  (30  to  50  mils) 
bonded  to  a  1-inch  thick  buffer  of  GE  type  151  fused  silica.  A 

* 

This  is  suggested  by  Royce  in  Reference  10. 
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reflecting  surface  was  put  on  the  fused  silica  by  vapor  depositing 
aluminum  on  the  front  surface  of  the  fused  silica  which  was  then 
bonded  to  the  aluminum.  The  interface  velocity  was  then 
monitored  with  the  laser  interferometer.  The  electron  beam 
machine  diagnostics  were  not  complete  enough  to  allow  a  detailed 
reduction  of  the  data.  However,  the  laser  interferometer  data 
are  interesting.  One  of  the  reduced  traces  is  shown  in  Figure 
33.  There  are  two  distinct  peaks  in  the  stress-time  profile. 

All  of  the  interferometer  data  taken  at  the  high  fluences  show 
evidence  of  the  second  peak  in  stress  time.  The  second  peak 
was  generally  of  larger  relative  amplitude  than  shown  in 
Figure  33  and  occurred  at  a  time  such  that  the  second  pulse 
could  not  be  explained  by  the  reverberation  of  the  initial 
stress  pulse  within  the  sample.  The  separation  in  time  of  the 
tv/o  peaks  was  approximately  the  same  as  predicted  by  the 
numerical  calculations.  While  the  experiments  are  not 
definitive,  they  do  provide  some  supporting  evidence  of  the 
existence  of  the  double  wave  structure. 

The  Chart  D  Analytic  Equation  of  State,  Reference  15  has 
not  been  included  in  the  comparisons.  This  EOS  model  is  a 
multi-constituent  model  which  may  provide  a  reasonable  alterna¬ 
tive  to  the  EOS  models  considered  in  this  report. 
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SECTION  6 
CONCLUSIONS 


The  objective  of  tnis  program  was  to  experimentally  and 
theoretically  investigate  the  radiation-induced  impulse  and 
stress  generation  characteristics  of  aluminum  for  cases  where  a 
significant  portion  of  the  blowoff  material  was  in  the  liquid- 
vapor  state. 

Experimental  data  were  obtained  in  the  laboratory  for  con¬ 
ditions  in  which  significant  amounts  of  the  blowoff  material 
were  in  the  liquid-vapor  state.  Impulse  measurements  were 
obtained  with  a  ballistic  pendulum  for  peak  doses  up  to 
3800  cal/gm.  Manganin  gauges  and  the  laser  interferometer  were 
used  to  obtain  stress  time  records  adjacent  to  the  region  of 
energy  deposition. 

The  PUFF  V,  GRAY,  and  PHILCO-FORD  equation-of-state  models 
were  compared  for  aluminum.  The  experimental  data  were  compared 
with  the  predictions  of  these  equation-of-state  models  plus  the 
LINEAR-FLUID  model.  The  PHILCO-FORD  equation-of-state  model 
gave  the  best  data  representation.  The  PUFF  V  model  and  the 
LINEAR  FLUID  model  were  found  to  give  approximately  the  same 
results  for  the  same  input  conditions. 

Examination  of  the  comparisons  of  the  data  with  the 
equation-of-state  models  provided  evidence  that  the  material 
response  of  aluminum  for  the  energy  loading  conditions  obtained 
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in  this  program,  in  which  a  significant  portion  of  the  blowoff 
was  in  the  liquid-vapor  state  at  relatively  low  pressures,  was 
not  dominated  by  mixed  phase  (liquid-vapor)  effects.  This  was 
because  the  maximum  pressure  of  the  liquid-vapor  state  is  only  a 
few  percent  of  the  maximum  pressure  due  to  energy  deposition. 

The  impulse  in  aluminum  for  the  experimental  conditions  obtained 
in  the  program  seems  to  be  dominated  by  the  behavior  of  the 
equation-of-state  of  aluminum  in  the  liquid  region  (see 
Figure  16  for  definition  of  liquid  region)  of  the  phase  plane. 

The  laser  interferometer  data  for  the  high-fluence  condi¬ 
tions  and  the  predictions  of  the  GRAY  equation-of-state  model 
suggest  the  possible  existence  of  a  double  wave  structure  in 
the  propagated  stress  pulse.  This  phenomenon  was  attributed  to 
the  recompression  of  the  liquid  layer  by  the  overlying  vapor 
layer  after  the  initial  drop  to  zero  pressure  in  the  liquid. 

With  the  installation  of  an  external-magnetic-field  beam- 
control  technique  on  the  pulsed  electron  beam  accelerators 
larger  areas  and  higher  peak  deposited  energies  have  become 
available.  The  larger  areas  will  lead  to  smaller  error  bars  on 
the  data  (Figure  26)  and  longer  one-dimensional  read  times. 
Specially  designed  experiments  could  be  used  to  investigate  the 
questions  which  remain  unanswered  after  this  program.  Some  of 
these  remaining  questions  are 

a.  Is  the  pressure-energy  coupling  coefficient  of 
reference  density  aluminum  at  internal  energy 
densities  above  2500  cal/gm  0.67  as  predicted  by  GRAY 
or  r-2.5  as  predicted  by  the  PHILCO-FORD  model? 

b.  Does  the  double  wave  structure  really  exist  for 
high  enough  fluences? 
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c.  Can  the  melting  and  vaporization  phase 
transformations  really  occur  on  the  tens  of 
nanosecond  time  scale  and  what  is  their  effect 
on  the  material  response? 

d.  What  is  the  equation  of  state  of  aluminum 
around  the  critical  region? 

This  program  has  concluded  that  the  material  response 
obtained  on  solid  aluminum  w;-.s  not  dominated  by  liquid-vapor 
effects.  However,  this  is  not  expected  to  be  the  case  for 
porous  metals .  In  the  case  of  porous  metals  the  material 
response  may  be  dominated  by  the  behavior  of  the  material  in 
the  liquid-vapor  regime.  The  importance  of  this  mixed  phase 
region  will  increase  as  the  porosity  increases.  Thus,  if  the 
material  response  of  porous  metals  is  to  be  understood,  a  well 
planned  theoretical  and  experimental  investigation  to  determine 
the  equation  of  state  of  metals  around  the  critical  region  must 
be  carried  out. 
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APPENDIX  A 

APPLICATION  OF  CALORIMETERS  AND  FARADAY  CUPS 
AS  FLUENCE  DIAGNOSTICS 


The  two  most  commonly  used  devices  to  diagnose  the  total 
energy  deposited  in  targets  by  electron  beams  are  the  total 
absorbing  carbon  calorimeter  and  the  Faraday  cup  charge  collector. 
The  topic  of  this  discussion  is  not  to  discuss  their  designs  but 
to  point  out  benefits  and  disadvantages  of  both  devices. 

Segmented  calorimeters  provide  the  experimenter  with  the 
most  reliable  existing  gauge  of  fluence  uniformity.  The  segments 
have  been  made  small  enough  to  fit  many  into  irradiation  areas 
which  are  typical  of  today's  electron  beams.  Assuming  that  all 
the  incident  beam  energy  is  deposited  in  the  calorimeter,  this 
device  becomes  a  useful  tool  in  measuring  the  incident  fluence 
and  its  spatial  distribution  provided  the  deposited  energy  stays 
in  the  calorimeter  long  enough  to  communicate  the  accompanying 
temperature  rise  to  the  attached  thermocouples.  Since  the 
thermocouples  must  be  attached  outside  the  energy  deposition 
region,  the  calorimeter  is  usually  slow  to  react.  Thus,  if  a 
fraction  of  the  deposited  energy  is  suddenly  removed  before  its 
temperature  rise  is  "sensed”  by  the  thermocouple,  the  resulting 
fluence  map  will  be  misleading. 

Open-shutter  photographs  of  attempted  high- fluence  measure- 
meants  using  segmented  graphite  calorimeters  (Figure  A-l)  as  well 
as  active  area  erosion  suggest  that  (hot)  material  is  removed 
during  fluence  measurements  for  which  the  graphite  sublimation 
energy  and/or  spall  strength  is  exceeded. 

Neglecting  removal  by  spallation  (which  may  be  aggravated 
by  moisture  absorbed  in  the  porous  graphite) ,  vapor  blowoff  from 
the  active  area  of  a  calorimeter  is  capable  of  carrying  away 
substantial  energy  which  cannot  be  communicated  to  the  thermo¬ 
couples  monitoring  temperature  changes  in  the  calorimeter  array. 
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Graphite  sublimation  energy  appears  to  be  approximately 
1500  cad/gm.  Whenever  graphite  calorimeters  are  used  to  measure 
fluences  in  excess  of  that  dose,  the  measured  fluence  is  likely 
to  be  less  than  the  deposited.  The  following  example  illustrates 
the  extent,  of  the  under  measurement  for  doses  in  excess  of 
1500  cal/gm. 


The  typical  normalized  peak  dose,  E 


max 


for  beams  with  mean 


electron  energy  of  260  keV  is  26  cal/gm/ (cal/cnr) .  With  a 
graphite  sublimation  energy,  Esub,  of  1500  cal/gm,  the  deposited 
fluence,  <J> 


Dep' 


fcr  which  graphite  just  begins  to  vaporize  is 


E  u 
sub 

.  max 


1500 
~7S~  = 


58  cal/cm 


2 


o 

Thus,  for  deposited  fluences  less  than  58  cal/gm  (at  260  keV) 
the  graphite  calorimeter  can  be  expected  to  measure  all  of  the 
deposited  fluence. 


Figure  A-2  illustrates  the  situation  for  deposited  fluences 
2 

above  58  cal/cm  . 


In  this  case,  graphite  will  be 
depth,  r,  at  which  the  dose  is  1500 


removed  by  vaporization  to  a 
cal/gm — i.e. , 


ESub 


<fc_  E  (r) 
Dep  o'  7 


The  area  under  the  non-normalized  deposition  profile  is  the 
deposited  fluence.  The  deposited  energy  corresponding  to  the 
unshaded  area  in  Figure  A-2  is  carried  away  by  the  vaporized 
material.  The  remaining  energy  (shaded  area)  is  communicated  to 
thermocouples  which  monitor  each  block  of  the  array  (Reference  1) . 
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From  the  geometry  in  Figure  A-2, 


ESub 

P(R-r) 


ib 

Pep 


so  that  the  apparent  fluence  measured  in  the  calorimeter  is 


vMeas 


1 

2 


P(R-r) Egufc 


PR(ESub)2  .  . 

^7^ for  ^ 


>58  cal/cm 


2 


In  general,  graphite  calorimeter  fluence  measurements  per¬ 
formed  with  beams  typified  by  a  260  keV  mean  electron  energy  will 
results  in  measured  fluence. 


2 

<f>  ,  for  <Jl,  <58  cal/cm 

Dep  Dep 


<P 


p*<Bsub> 


Me  as 


2  6  E 
vDep  o 


,  for  a  >58  cal/cm2 
Max  rDep 


The  measured  fluence  is  plotted  against  deposited  fluence 
in  Figure  A-3  for  mean  electron  energies  of  220,  260,  and  300  keV. 
The  maximum  in  each  of  the  three  curves  identifies  the  critical 
deposited  fluence  beyond  which  calorimeter  graphite  is  removed 
by  vaporization.  Thereafter,  measured  fluence  decreases  with 
increasing  deposited  fluence. 
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As  mentioned  earlier,  the  1500  cal/gm  limit  is  an  upper 
bound  beyond  which  graphite  calorimetry  should  not  be  used. 
Front-surface  spall  and  the  eruptive  effect  when  moisture  within 
the  graphite  pores  changes  phase  suggest  that  graphite  calori¬ 
meters  may  not  be  accurate  above  3000  cal/gm. 

Although  this  limitation  is  disappointing,  the  graphite 
calorimeter  is  none  the  less  a  useful  tool  for  low-dose  work. 

It  is  still  the  only  diagnostic  capable  of  quantitative  spatial 
resolution  of  the  deposited  energy.  Furthermore,  since  calori¬ 
metry  data  is  fast  to  reduce,  it  provides  on-the-spot  fluence 
data  during  the  course  of  an  experiment. 

Measurements  of  the  total  beam  energy  under  high-dose  con¬ 
ditions  are  performed  using  the  outputs  of  the  diode  voltage 
monitor  and  a  Faraday  cup  as  detailed  in  Reference  1.  Integra¬ 
tion  under  the  power  curve  yields  the  to^al  beam  calories  de¬ 
posited  in  the  Faraday  cup 


H  -  37186  / V(t>  dt 


Because  the  Faraday  cup  measures  current,  it  is  insensitive 
to  vaporization  or  front-surface  spall  of  its  collector.  The 
penalty,  in  this  case,  is  the  lengthy  and  tedious  process  of 
converting  oscilloscope  traces  into  power  curve  integrals. 

The  disadvantage  of  both  the  carbon  calorimeter  and  the 
Faraday  cup  techniques  is  that  neither  accounts  for  the  energy 
back-scattered  out  of  the  measurement  device.  Tne  remainder  of 
this  section  will  be  devoted  to  methods  of  accounting  for  the 
back-scattered  energy,  at  least,  to  a  reasonable  approximation. 


Ill 


It  will  be  argued  later  that  the  waveform  (shape)  of  the 
incident  current  is  similar  to  the  current  absorbed  by  the 
Faraday  cup.  This  circumstance  makes  the  input  to  the  Monte 
Carlo  electron  deposition  code  not  an  unreasonable  one  even  in 
view  of  backscatter.  The  input  requires  only  the  current  time 
dependence,  not  its  absolute  magnitude,  in  addition  to  the 
accelerating  voltage  waveform  and  the  mean  incidence  angle  of 
the  beam  electrons.  In  short,  the  code  "sees"  a  primary  current 
waveform  (from  the  Faraday  cup)  not  unlike  the  (unmeasured) 
incident  primary  current,  together  with  the  voltage  history  which 
accelerated  the  (unmeasured)  incident  primary  current.  In 
addition  to  its  main  function  of  calculating  the  energy  depo¬ 
sition  profile,  the  code  calculates  the  fraction  of  incident 
energy  back-scattered,  <£>  ,  as  well  as  the  percent  of  incident 
electrons  back-scattered,  <a>T.  Both  <ct>x  and  <3>T  are  given 
as  averages  over  the  entire  beam  pulse.  Figure  A-4  demonstrates 
that  both  quantities  have  weak  dependence  on  mean  electron 
energy  but  strong  dependence  on  the  mean  incidence  angle  of  the 
beam  electrons  in  the  range  of  interest  for  this  program. 

Since  calorimeters  measure  the  absorbed  energy  over  the 
entire  pulse,  the  pulge-averaged  <B>  is  appropriate  and  the 
back-scatter  correction  is  trivial  once  the  mean  incidence 
angle  is  known.  If  H.  is  the  energy  incident  or.  a  calorimeter, 
Hcal'  the  absorbed  and  Href  the  back-scattered  energy,  then 


H.  =  H  _  +  H 
me  cal  ref 


H 


<B> 


ref 


T  “  H, 


me 


so  that 


H. 


me 


H 

cal 


Knowledge  of  the  mean  incidence  angle  defines  (see 

Figure  A-4)  and  the  back-scatter  correction  is  completed. 

Since  the  voltage  time  history  of  the  back-scattered 
electrons  is  unknown,  a  <0>  -correction  cannot  be  applied  to 
the  power  curve  integral.  However,  an  attempt  can  be  made  to 
account  for  backscatter  using  the  fraction  of  electrons  back- 
scattered. 

V(t)  in  the  power  curve  integral  needs  no  correction  since 
it  already  describes  the  accelerating  voltage  of  the  (unmeasured) 
incident  primary  current.  The  inadequacy  arises  when  the  primary 
current  absorbed  in  the  Faraday  cup,  IF(,(t)  •  is  used  in  lieu 
of  the  unmeasured  incident  primary  current,  If  iref^ 

represents  the  current  associated  with  the  back-scattered  electrons 
then 


W*'  "  +  Iro£<t> 


where  each  of  the  three  currents  could  have  a  different  time 
dependence.  If  a(t)  is  the  fraction  of  incident  electrons 


reflected  at  any  time,  t,  in  the  beam  pulse,  then 


a(t)  = 


Wt)  dt 


hncM  "  IFC(t>  +  “(t)  hncM 


where 


^nc'4' 


IpC(t) 

l-o(t) 


This  general  statement  is  still  not  useful  since  only  IpC(t) 
is  known  from  the  Faraday  cup  measurement. 

If  it  can  be  shown  that  a(t)  is  reasonably  constant  over 
the  beam  pulse,  the  backscatter  correction  is  made  possible  by 
using 


a(t)  =  constant  =  <a> 


4.186  (l-<a>  ) 


y  v< t)ipc(t)  dt 


A  constant  a(t)  requires  that  the  three  currents  I^nc(t), 
IpC(t),  and  all  have  the  same  waveform  and  differ  only 

in  amplitude.  The  similarity  of  time  dependence  (shape)  cannot 
be  directly  verified  since  l£nc (t)  is  not  measured  at  the  Faraday 
cup  location.  However,  the  unmeasured  incident  primary  current, 
I-  (t),  is  generated  in  the  diode  where  the  diode  current  is 

me 

monitored  by  a  B-probe  (Reference  1) .  Figure  A-5  illustrates 

that  the  current  absorbed  in  the  Faraday  cup,  IpC(t),  has  a 

shape  similar  to  that  of  the  diode  currents  which  is  the  initial 

shape  of  the  incident  primary  current,  i^c^t).  Since  even  the 

oscillatory  structure  is  the  same  in  the  two  waveforms,  it  is 

reasonable  to  assume  that  (to  a  good  approximation)  the  diode 

current  changes  only  in  amplitude  and  beam  front  erosion  as  it 

propagates  toward  the  Faraday  cup  to  become  the  unmeasured 

I.  (t) .  Thus,  a(t)  appears  to  be  very  nearly  constant  over  the 
me 

beam  pulse  and 


H  • 
me 


f v(t)iFC(t)  dt 


becomes  a  reasonable  backseat ter  correction  provided  the  mean 
incidence  angle  of  the  beam  electrons,  <8>,  is  known  in  order  to 
define  the  appropriate  <ct>T  in  Figure  A- 4. 

Reference  3  gives  an  expression  of  <8>  in  terms  of  the 
primary  beam  current,  I  ,  its  relativistic  By  and  the  current 
neutralization  factor,  (1-fm) ,  defined  as  the  ratio  of  net  to 
primary  current  amplitudes.  The  equation,  derived  from  magnetic 


current,  relative  units 


Time,  nsec 


Figure  A- 5  Comparison  of  current  waveforms. 


pressure  balance,  is 


The  backscatter  correction  given  above  complicates  matters 
since  now 

■^FC 

T  =  T  =  — _ — 

PR  inc  l-«x>T 

2 

whereby  the  tan  <0>  expression  becomes  one  equation  m  two 
unknowns — i.e.,  <0>  is  needed  to  define  <a>T  and  <a>T  is  needed 
to  calculate  <0>. 

The  only  recourse  is  to  adopt  an  iterative  procedure  taking 
advantage  of  the  fact  that  tan2  <0>  is  a  slow-varying  function 
of  <a>x. 

The  backscatter  corrections  described  in  this  appendix  were 
applied  to  the  data  taken  as  part  of  this  program. 
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SECTION  B-l 

LASER  VELOCI  METER 


The  Physics  International  laser  interferometer  system  is 
designed  to  measure  the  response  of  materials  to  electron  beam 
deposition  or  mechanical  shock  loading.  Basically,  it  records 
the  time-velocity  and  time-position  profile  of  the  rear  surface 
of  a  sample  under  test.  This  surface  may  be  either  free  or  in 
contact  with  a  suitable  buffer  material  such  as  fused  silica  or 
lucite.  To  date,  measurements  have  been  made  on  aluminum,  sparesyl, 
fused  silica,  Comrad  I  fiber  bundles,  and  several  phenolics  with 
excellent  results.  Electron  beam  conditions  range  from  front 
surface  deposition  to  uniform  deposition  throughout  the  thick¬ 
ness  of  the  sample.  Measured  velocities  have  ranged  between 
5  m/sec  and  650  m/sec  with  accuracies  on  the  order  of  0.1  percent. 

An  outline  of  the  optical  configuration  used  in  this  system 
is  shown  in  Figure  B-l.  The  beam  from  an  argon  ion  laser  is  re¬ 
flected  from  the  rear  surface  of  the  sample.  As  the  surface 
accelerates  under  shock  loading  conditions,  the  reflected  beam 
is  given  a  phase  shift  and  a  Doppler  shift  due  to  the  surface 
displacement  and  velocity,  respectively.  Optical  interference 
is  used  to  measure  the  size  of  these  shifts.  Part  of  the  return 
beam  is  combined  with  a  static  reference  beam  for  displacement 
information  and  the  remainder  of  the  beam  is  split;  part  of  it 
is  delayed  prior  to  recombination  to  measure  velocity.  Upon 
recombination,  the  relative  intensity  of  each  beam  is  propor¬ 
tional  to  the  cosine  of  the  relative  phase  angle.  The  intensities 
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Figure  B-l  Tandem  mode  interferometer  system  offering  simultaneous 
displacement,  velocity,  and  acceleration  information. 

are  monitored  by  two  photomultiplier  tubes  and  data  appear  as  a 
series  of  fringes  (360  degrees  phase  shifts). 

In  .the  displacement  mode,  each  fringe  indicates  a  change 
in  the  sample's  rear  surface  position  of  2.44  x  10**5  centimeters 
along  an  axis  normal  to  that  surface.  In  the  velocity  mode,  each 
fringe  corresponds  to  a  change  in  surface  velocity  of  18  m/sec. 

This  correspondence  between  fringes  and  change  in  velocity  in  the 
latter  mode  can  be  altered  for  any  given  experiment  depending 
upon  the  anticipated  rear  surface  behavior.  By  changing  the 
length  of  the  delay  leg  the  correspondence  may  be  varied  from 
10  m/sec/fringe  to  over  1000  m/sec/fringe. 
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1.1  TECHNICAL  DISCUSSION 


The  laser  system  in  Figure  B-l  can  be  broken  down  into  three 
basic  components— transmission- reception  unit,  displacement 
interferometer,  and  velocity  interferometer.  The  transmission- 
reception  unit  (Figure  B-2)  consists  of  the  laser,  a  spatial 
filter,  a  variable  beam  splitter,  a  polarizing  beam  splitter, 
and  a  quarter-wave  plate.  A  1.2  watt  light  output  is  supplied 
by  the  laser  at  a  wavelength  of  4880  A.  To  increase  the  signal- 
to-noise  ratio  the  beam  is  first  put  through  a  spatial  filter 


Expanded  view  OS  spatial 

filter  (beam  diameter  exaggerated) 
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Reference  To 

interferometers 


Spatial  filter  beam  Polariting  Quarter 

aditter  boam  splitter  vave  plate 


Figure  B-2  Transmission-reception  unit 


consisting  of  two  biconvex  lenses  with  a  pinhole  stop  placed 
between  them  (Figure  B-2) .  Since  the  output  of  the  laser  is  mono¬ 
chromatic  and  collimated,  the  beam  is  focused  to  a  point  by  the 
first  lens.  The  pinhole  stop  is  placed  at  the  focal  point  of 
the  lens  where  the  beam  passes  through  to  be  re collimated  by  the 
second  lens.  Noise  light  generated  in  the  laser  tube  is  not  well 
collimated,  does  not  become  focused,  and  is  blocked  by  the  pinhole 
stop.  The  focal  lengths  of  the  two  lenses  are  in  the  approximate 
ratio  of  2:1  which  affects  a  similar  reduction  in  beam  diameter. 
This  diameter  reduction  was  deemed  necessary  due  to  the  long 
path  lengths  over  which  the  beam  must  be  transmitted  in  the 
experimental  configuration.  An  increase  in  signal- to-noise 
ratio  of  approximately  3:1  was  noted  due  to  the  spatial  filter. 
Alignment  of  the  detection  system  was  also  simplified  since  it 
was  no  longer  necessary  to  search  for  a  non-noisy  spatial  position 
of  the  return  beam  to  interrogate  the  photomultiplier  tube. 

From  the  spatial  filter,  the  beam  passes  through  a  variable 
beam  splitter,  a  polarizing  beam  splitter,  a  quarter  wave  plate, 
and  is  then  transmitted  to  the  sample.  The  variable  beam 
splitter  diverts  a  small  portion  of  the  light  energy  to  be  used 
as  a  reference  beam  for  displacement  mode  information.  The  output 
of  the  laser  is  polarized  and,  due  to  the  orientation  of  its 
plane  of  polarization  relative  to  the  polarizing  beam  splitter, 
all  of  the  light  energy  is  transmitted  at  this  point. 

In  the  electron  beam  chamber  the  laser  beam  is  focussed 
to  a  10  micron  spot  on  the  sample's  surface  with  a  long  focal 
length  lens.  Since  the  surface  is  at  the  focal  point  of  the 
lens,  the  reflected  laser  beam  returns  along  its  incident  path. 
This  technique  simplifies  alignment  and  reduces  the  sensitivity 
to  sample  tilt.  Slight  tilt  on  the  part  of  the  sample  will  result 
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in  only  a  small  lateral  displacement  of  the  return  beam.  On  re¬ 
turning  to  the  polarizing  beam  splitter,  the  beam  is  reflected 
instead  of  being  transmitted  since  its  plane  of  polarization  has 
been  rotated  90  degrees  by  a  double  transit  through  the  quarter 
wave  plate.  Thus,  the  return  beam  is  prevented  from  re-entering 
the  laser  cavity  to  cause  instability.  From  here  the  return  beam 
is  passed  through  the  two  interferometers  and  its  information 
content  is  converted  to  the  outputs  of  the  two  photomultiplier 
tubes . 

Displacement  information  is  obtained  by  diverting  part  of 
the  energy  of  the  return  beam  with  a  second  variable  beam  splitter 
and  combining  it  with  the  reference  beam  as  shown  in  Figure  B-l. 
The  path  length  of  the  reference  beam  is  adjusted  with  a  double 
mirror  bank  to  be  approximately  equal  to  twice  the  distance  from 
the  laser  to  the  sample  under  inspection;  this  ensures  that  the 
path  length  remains  within  the  coherence  length  of  the  system. 
Since  the  laser  has  a  coherence  length  of  10  meters,  this  require¬ 
ment  is  not  singularly  critical.  From  the  displacement  interfero¬ 
meter,  the  beam  passes  through  a  diverging  lens  and  into  the 
photomultiplier  tube.  The  diverging  lens  allows  only  one  fringe 
to  be  interrogated,  eliminating  spatial  integration  effects. 

Also,  it  allows  the  input  light  intensity  to  be  adjusted  for 
maximum  signal  and  signal-to-noise  ratio.  A  typical  displacement 
mode  trace  along  with  the  velocity-time  profile  is  shown  in 
Figure  B-3. 

The  portion  of  the  return  beam  not  used  in  the  displacement 
mode  is  interrogated  by  the  velocity  interferometer.  Here  the 
Doppler  shift  is  detected  by  comparing  the  phase  of  the  light  at 
one  instant  of  time  with  phase  at  an  earlier  instant.  The  rela¬ 
tive  phase  is  then  directly  proportional  to  the  rear  surface 
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velocity.  A  velocity  interferometer  (Figure  B-4)  is  used  to 
affect  this  phase  comparison.  The  return  beam  is  split  at  the 
first  beam  splitter  and  half  of  the  light  is  directed  around  a 
delay  leg  to  be  recombined  with  the  direct  beam  at  the  second 
splitter.  A  double  concave  lens  diverges  the  recombined  beam 
that  is  then  detected  by  the  photomultiplier  tube.  As  the  rear 
surface  velocity  changes,  so  does  the  relative  phase  of  the  re¬ 
combined  beam.  For  a  shift  of  one  full  fringe,  a  minimum  change 
in  surface  velocity  of  10  m/sec  is  required.  This  relationship 
may  be  altered  where  very  high  rear  surface  velocities  and 
accelerations  are  anticipated  by  shortening  the  length  of  the 
delay  leg.  The  relationship  between  the  number  of  fringes  observed 
and  the  rear  surface  velocity  is  given  by  the  following  formula: 


where  u  is  the  free-surface  velocity,  A  is  the  wavelength  of  the 
laser  light,  C  is  the  speed  of  light,  D  is  the  length  of  the 
delay  leg,  and  N  is  the  number  of  fringes  observed  up  to  that 
point  in  time.  Figure  B-5  shows  a  typical  example  of  the  photo¬ 
multiplier  tube  output  and  associated  data  reduction  for  a  test 
on  reference  Sparesyl. 
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126 


SECTION  B -2 

OPTICAL  DETECTORS 


Photomultiplier  tubes  have  been  used  as  a  light  sensor  to 
detect  destructive  and  constructive  fringes  of  the  interfering 
laser  beams  and  to  convert  these  changes  in  light  intensity  to 
electrical  signals.  Three  RCA  PM  tubes  have  been  tested  in  the 
course  of  developing  the  Physics  International  laser  interfero¬ 
meter  into  an  effective  diagnostic  tool  for  the  study  of  dynamic 
behavior  in  materials.  They  are  models  1P28,  6810A,  and  C-70045C. 
Due  to  its  high  gain,  tube  model  6810A  was  found  to  be  overly 
sensitive  and  hard  to  stabilize.  For  this  reason,  the  other  two 
tubes  were  chiefly  used.  The  1P28A  model  was  used  for  bench 
testing  and  experiments  requiring  frequency  response  under  100  MHz; 
model  C-70045C  was  used  for  light  signals  of  higher  than  100  MHz 
frequency.  While  Model  1P28  was  used  regularly  in  both  CW  and 
gated-modes,  the  C-70045C  tube  was  mainly  incorporated  in  a 
pulsed  scheme  with  optional  CW  operation.  A  brief  description 
of  the  system  employing  the  fast  RCA  photomultiplier  tube  is 
given  below. 

B.2.1  PM  TUBE  AND  ELECTRICAL  CIRCUIT 

The  C-70045C  tube  is  an  experimental  PM  tube  design  from 
RCA.  It  is  a  side-on  device  having  an  S-20  response  (Figure  B-6). 
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Figure  B-6  Typical  spectral  response  characteristics 
of  RCA-C70045  photomultiplier  tube. 


Fourteen  dynodes  alternate  with  14  accelerating  electrodes  and 
a  50-ohm  coaxial  output.  The  factory  objective  sheet  gives  a 
risetime  of  0.5  nsec  and  a  gain  of  10^.  Displaying  the  anode 
response  of  an  occulted  tube  to  cosmic  ray  pulses  with  a 
Tektronix  7074  oscilloscope  (7A16  plug-in)  has  reflected  the  rise- 
time  limitation  of  the  scope  (T  =2.4  nsec)  rather  than  that 
of  the  tube  (Figure  B-7) .  The  frequency  response  of  the  C-70045C 
in  laser  interferometry  appears  to  be  limited  by  the  pulse  fall¬ 
time  which  was  found  to  be  1.5  to  2.0  nsec. 
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Note:  The  risetime  is  typical  of  the  oscilloscope  used 
(Tektronix  7704). 

Figure  B-7  Anode  output  pulse  of  an  occulted  RCA  C-70045C 
PM  tube  (V  =  5000  V) 

In  developing  electrical  circuits  for  the  PM  tube,  the 
findings  of  similar  experimental  programs*  have  been  studied 
and  adapted  where  they  were  found  applicable.  Consequently, 
several  modifications  to  the  manufacturer's  suggested  circuit 
have  been  made  to  avoid  overload  or  "run-away"  of  the  tube  at 
dc  light  levels  expected  to  be  attained  within  the  interferometer. 
The  most  important  modification  was  designed  of  a  gating  circuit 
to  operate  the  PM  tube  under  its  optimum  conditions ,  that  is , 
to  detect  low-level  light  pulses  of  short  duration.  The  circuit 
turns  on  the  tube  by  applying  a  voltage  pulse  between  the  photo¬ 
cathode  and  dynode  number  3,  both  floating  normally  at  the  same 
potential  in  a  turned-off  condition.  The  circuit  is  a  simple 
SCR  design  (Figure  B-8)  which  discharges  a  0.05  nF  capacitor, 
thereby  adding  1200  to  1600  volts  for  a  period  of  25  ysec  to  the 
aforementioned  divider  section. 

* - 

Sandia  Corporation  Report  SC-DR-66-2689 . 


130 


Minibox 


1.  Gate  output:  ~  -  U/5  V,*  for  PMT  cathode 

2.  +  V2:  typically  2000  volts  for  RCA.  C-70045C 

3.  Trigger  pulse:  gate  pulse  from  Tektronix  555 

Figure  B-8  Gating  circuit  used  with  PM  tubes  on  the 
laser  interferometer. 

The  gating  and  voltage  divider-circuits  are  designed  to 
operate  each  of  the  32  elements  at  the  factory  recommended  voltage 
The  voltage  divider  circuit  (Figure  B-9)  is  a  4  mA  resistor 
divider  with  interdynode  capacitors  mounted  on  the  socket.  The 
large  capacitors  are  used  to  maintain  constant  voltage  and  gain 
per  stage  when  the  tube  is  gated-on  in  the  presence  of  high  light 
levels  while  the  smaller  capacitors  are  necessary  to  stabilize 
the  gain  during  high  frequency  light  transients. 
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The  maximum  signal  currents  obtained  while  gating  without 
overloading  the  tube  were  70  to  80  mA  through  the  50-ohm  load 
resistor.  It  was  decided  that  biasing  the  anode  of  the  PM  tube 
(thus  further  increasing  the  gain  at  the  expense  of  possible 
decreased  frequency  response)  would  not  be  used. 

B.2.1  EXPERIMENTAL  SETUP 

A  typical  arrangement  for  displaying  interferometer  detector 
signals  is  shown  in  Figure  B-10.  The  photomultiplier  tube  with 
a  cylindrical  light  channel  incorporating  a  1-inch-diameter 
aperture  is  facing  the  laser  beam  and  is  turned  on  by  the  gating 
circuit.  Via  a  555  oscilloscope,  the  Pulserad  control  panel 
triggers  the  tube  7  psec  in  advance  of  the  electron  beam  pulse. 
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Figure  B-10  Block  diagram  of  a  typical  experimental  setup. 
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The  photomultiplier  anode  circuit  was  completed  through  an 
external  50-ohm  load  resistor;  due  to  the  signal  current  the 
voltage  drop  across  the  load  was  displayed  by  several  454  type 
Tektronix  oscilloscopes.  In  order  to  show  the  temporal  history 
of  the  signal,  the  sweeps  of  the  scopes  were  cascaded  by  using 
the  delayed  gate  output  of  a  scope  to  trigger  the  succeeding  one. 
The  first  oscilloscope  in  the  chain  was  triggered  typically  by 
the  B  probe  signal  of  the  Pulserad  after  selection  of  sweep 
speeds.  A  type  184  time  mark  generator  was  used  to  synchronize 
the  participating  oscilloscopes.  Normally,  a  two  division  over¬ 
lap  was  maintained  between  two  displays  to  utilize  only  the  most 
linear,  central,  8-centimeter  section  of  each  oscillograph. 

During  the  pre-shot  alignment  of  the  interferometer,  the 
photomultiplier  tube  was  repetitively  gated  (60  Hz)  to  continuously 
view  the  laser  beam  modulated  by  an  oscillating  mirror  of  the 
system.  For  this  quasi  CW  test,  4,000  volts  were  maintained 
between  dynode  number  3  and  the  anode,  with  1,200  volts  switched 
onto  the  divider  section,  thus  turning  on  the  photocathode. 

Light  intensity  was  adjusted  to  slightly  overload  the  PM  tube 
with  maximum  ac  to  dc  signal  ratio.  Finally,  the  detector  system 
was  readied  for  the  data  shot  by  switching  to  the  single  gate 
scheme  with  the  expected  voltage  pulse  between  dynode  number  3 
and  the  photocathode  increased  from  1,200  to  1,600  volts.  Typical 
signal  levels  of  2  to  3  volts  were  obtained  with  signal-to-noise 
ratios  in  excess  of  20:1  using  this  system. 
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SECTION  B-3 

RESULTS  AND  CONCLUSIONS 


The  advantages  of  the  laser  interferometer  as  a  diagnostics 
tool  are  several.  An  optical  probe  applies  no  loading  to  the 
surface  under  test  for  free  rear  surface  measurements.  In 
contrast  to  quartz  gauges /  there  is  no  averaging  of  the  stress 
over  an  area  because  of  the  small  spot  size  used  (approximately 
10  n  diameter) .  There  is  also  no  intrinsic  readtime  limitation 
as  with  quartz  gauges;  (time  windows  of  up  to  25  ptsec  have  been 
used) .  Alignment  of  the  system  has  been  simplified  and  turn¬ 
around  time  of  as  low  as  20  minutes  between  shots  has  been 
achieved.  Presently  the  laser  interferometer  is  being  used  in 
conjunction  with  a  ballistic  pendulum  for  making  simultaneous 
measurements  of  impulse  and  stress. 


Etched  manganin-foil  gauges  {3/4-mil  thickness)  and 
manganin  wire  gauges  have  been  used  in  the  past  at  PI.  The  foil 
gauges  (1.5-ohm  impedance)  are  proposed  for  the  tests  since  they 
offer  the  advantage  of  being  sandwiched  between  two  pieces  of 
test  material  and  are  readily  available  in  large  quantity.  Use 
of  the  wire  gauges  (about  7-ohm  impedance)  potted  in  C-7  epoxy 
necessitates  a  correction  for  differences  in  shock  impedance 
between  the  epoxy  and  the  sample  material.  If  the  shock  imped¬ 
ances  of  any  test  materials  are  well  known,  this  correction  is, 
of  course,  not  a  problem.  Signal-to-noise  ratios  are  suf¬ 
ficiently  large  to  enable  recording  of  a  1-kbar  stress  pulse 
in  the  foil  gauges  (shown  in  the  following  section) .  In  prac¬ 
tice  we  have  seen  very  little  difference  in  signal-to-noise 
ratio  between  the  foil  and  wire  gauges. 

C.l  MANGANIN  GAUGE  PRESSURE  RECORDING  SYSTEM— EXPERIMENTAL 
SETUP 

Figure  c-1  is  a  schematic  diagram  of  the  manganin  recording 
system,  similar  to  that  used  by  Stanford  Research  Institute. 

The  power  supply,  which  provides  a  nearly  constant  current  of 
1  to  5  amps  to  the  gauge,  is  triggered  approximately  70  usee 
prior  to  beam  firing  by  the  delay  trigger  generator.  The  "turn 
on"  voltage,  V,  produced  across  the  active  gauge  resistance, 

R,  is  monitored  on  one  trace  (channel  A)  of  a  Tektronix  555 
oscilloscope  (Figure  C-2) .  Arrival  of  a  plane  shock  wave  at 
the  active  gauge  element  causes  a  change  in  its  resistance, 

AR,  which  results  in  a  change  in  the  voltage,  AV,  monitored  on 
channel  A.  The  relative  voltage  change,  AV/V,  is  typically  <1 
percent  under  lowrstress  conditions  and  precludes  accurate 
measurement  of  AV  on  the  channel  A  oscilloscope  trace.  Record¬ 
ing  of  AV  is  accomplished  by  subtraction  of  the  turn-on  voltage 
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D  Probe  (Diode  current) 
output  attenuated  10:1 


^  Trigger  for  differential 
amplifiers  (Type  W  or  1A5) 
recording  channel  A  minus 
channel  B 


Figure  C-l  Schematic  of  manganin  measuring  svstem. 
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(V)  from  the  gauge  output  (v  +  AV)  using  a  differential  ampli¬ 
fier  {a  Type  W  plug-in  unit)  by  splitting  the  gauge  output  into 
two  channels,  with  channel  A  going  directly  to  the  scope.  The 
difference  A-B  is  monitored  on  the  other  trace  of  the  555  oscil¬ 
loscope,  and  this  trace  is  triggered  at  the  time  of  beam  firing. 
Figure  C-3  shows  the  time  history  of  this  A-B  signal  and  indi¬ 
cates  the  portion  displayed  on  the  oscilloscope.  The  term  t 
is  the  time  at  which  the  power  supply  is  turned  on.  The  term 
t^  corresponds  to  the  time  for  the  initial  signal  to  travel 

through  the  delay  line  =  9.2  usee  for  6000  ft  of  RG  213/U 

cable) .  The  term  t2  is  the  time  it  takes  for  the  system  to 

stabilize  (t2  =  60  usee  for  6000  ft  of  cable) ;  t^  is  the  time 

of  the  signal  arrival  of  channel  A  at  the  scope.  This  recording 
system  reads  AV(t)  until  time  t4  =  t3  +  9.2  usee  for  a  6000  ft 
delay  cable.  After  t^  the  signal  becomes  s(t)  =  AV(t)  -  AV 
(t  =  9.2  usee).  Note  that  the  dc  resistances  of  channels  A  and 
B  must  be  equal  (a  variable  resistor  was  used,  R5,  to  balance 
out  the  dc  resistance  of  the  delay  cable) .  With  the  present 
system,  V  is  measured  on  one  trace  and  AV  on  another.  A  dc 
circuit  analysis  gives  the  relationship  between  AV/V  and  AR/R 
(Figure  C-4)  from  which  the  stress  can  be  calculated. 

In  order  to  operate  the  manganin  gauge  in  the  electron-beam 
environment,  it  must  be  electrically  shielded.  This  was  done 
by  placing  the  gauge,  target,  and  cables  inside  a  4-inch-id 
variable-length  copper  cylinder  with  1/16-inch-thick  walls  that 
bolts  to  the  aluminum  drift-chamber  door  (Figure  C-5).  The 
front  face  is  a  1/2-inch-thick  aluminum  plate  with  an  orifice 
just  large  enough  to  admit  the  end  of  the  guide  cone.  The 
target  is  held  against  the  orifice  by  a  brace.  In  the  case  of 
nonmetallic  targets  it  is  necessary  to  place  a  grounded  piece 
of  0.001-inch-thick  aluminum  foil  between  the  target  and  the 
gauge  to  shield  the  gauge  element  from  noise. 
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AV/V  at  the  scope  versus  Af./R  at  the  gauge 


Drift  chamber 


Current  in 


Current  out 


Current  out 


Figure  C-5  Detai]  of  shielding  configuration 


2. 


Results 


Results  of  manyanin  tests  in  a  two-dimensional  weave 

material  are  presented  here  to  exhibit  signal-to-noise  ratio 

versus  time,  recording  time,  and  gauge  sensitivity  under  actual 

operating  conditions.  Some  1.5-ohm  foil  gauges  were  used  with  a 

1-mil  aluminum  shielding  foil  placed  between  the  front  sample 

section  and  the  gauge  element.  One-half-mil  Mylar  was  used  as 

electrical  insulation  on  each  side  of  the  gauge.  The  gauge  was 

sandwiched  between  a  1/4-inch  piece  of  the  test  material  (v/hich 

2 

was  irradiated  with  250  keV  electrons  at  about  45  cal/cm  uniform 
2 

over  2.5  cm  )  and  a  1/2-inch- thick  piece  of  test  material.  The 
shielding  canister  was  also  used. 


Figure  C-6  shows  the  recorded  gauge  output.  Noise  level  is 
seen  to  decay  to  about  10  mV  in  0.5  ysec,  and  drops  to  <  5  mV 
after  1.5  ysec.  Noise  level  remains  at  this  level  until  9  ysec 
at  which  time  noise  arrives  through  to  delay  line,  giving  an 
effective  reading  time  of  about  8  nsec. 

The  stress  pulses  shown  in  Figure  C-6  are  from  the  same 
shot.  Figure  C-6 a  shows  the  shock  pulse,  on  scale,  followed  by 
a  low-stress,  long-duration  pulse  which  is  more  clearly  visible 
in  Figure  C-6b  (higher  scope  sensitivity  scale) .  Peak  pressure 
of  the  shock  pulse  is  approximately  8  kbar,  while  the  following 
pulse  peaks  at  approximately  2  kbar.  Signal-to-noise  ratios 
for  the  whole  stress  record  are  seen  to  be  about  10  to  1. 

Integration  of  pressure  versus  time  gives  an  impulse  of 
5.0  ±  1.0  ktap  (the  uncertainty  arises  from  a  possible  baseline 
drift  at  late  times) .  Independent  measurement  of  total  impulse 
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under  conditions  similar  to  the  shot  illustrated  in  Figure  C-6 
show  an  impulse  of  4.2  ±  0.3  ktap.  The  conclusion  of  this  work 
was  that  impulse  was  delivered  to  the  material  by  the  initial 
shock  pulse  (short  duration)  and  by  a  long-duration  pulse 
caused  by  front- surface  blowoff. 

Other  examples  of  manganin  results  are  shown  in  Figures 
C-7  and  08.  These  data  were  collected, as  part  of  a  program 
for  AVCO,  with  manganin  wire  gauges.  The  spike  visible  in 
Figure  C-7  at  about  9  psec  is  the  arrival  of  noise  through  the 
delay  cable.  As  explained  in  the  preceding  section  the  pressure 
time  profile  can  oe  reconstructed  from  the  gauge  output  for 
times  longer  than  9  ysec.  In  this  particular  material  no 
double-pulse  structure  was  evident. 

3.  Simultaneous  Impulse  and  Stress-Time  Measurements 

The  manganin-gauge  technique  described  above  was  combined 
with  the  ballistic  pendulum  measurement  of  impulse.  The  gauge 
was  positioned  in  the  pendulum  bob.  The  gauge  power  and  signal 
leads  were  run  out  of  the  bob,  through  the  pendulum  rod  and  axle 
to  a  terminal  board  on  the  axle  (Figure  C-9).  Damping  of  the 
pendulum  motion  was  avoided  by  continuing  the  connections  from 
there  by  way  of  loosely  draped,  thin  wires  to  a  stationary 
terminal  board  on  the  pendulum  support.  RG  58  cables  were  used 
to  complete  the  circuit. 

The  gauge  leads  were  shielded  from  the  RF  environment  in 
the  drift  chamber.  The  pendulum  bob,  rod  and  axle  provided 
shielding  up  to  the  decoupler  (Figure  C-10).  The  decoupler  and 
the  cables  up  to  the  beam  chamber  wall  were  shielded  by  a 
flexible,  copper  pipe  (see  foreground  of  Figure  C-10) . 
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When  the  aluminum  samples  are  examined  after  irradiation  by 
the  pulsed  electron  beam,  the  characteristic  appearance  of  the 
craters  is  observed  tc  change  rather  dramatically  as  a  function 
of  the  mean  electron  energy  and  beam  intensity.  For  low  (200  to 
300  keV)  electron  energy  beams,  the  craters  at  low  fluence  are 
clean  and  the  floor  exhibits  a  granular  texture  similar  in  appear¬ 
ance  to  rear  surface  spalls,  but  somewhat  smoother  (Figure  D-l 
(a) ) .  As  the  fluence  is  increased,  the  crater  depth  increases 
and  droplets  of  solidified  melt  are  observed  on  the  crater  floor 
(Figure  D-l  (b) )  .  As  the  fluence  is  increased  further,  more 
and  more  solidified  melt  is  found  in  the  crater,  until,  at  the 
highest  fluence,  the  crater  is  obscured  by  the  solidified  melt 
which  exhibits  evidence  of  extensive  radial  flow  (Figure  D-l 
(c) ) .  Using  higher  energy  beams  (>  600  keV)  the  craters  are 
deep  and  clean,  even  at  high  fluences. 

The  observed  variation  in  crater  appearance  and  the  extent 
of  mass  removal  can  be  qualitatively  explained  by  considering 
the  basic  physical  processes  that  produce  front-surface  mass 
loss.  To  simplify  the  following  discussion,  effects  due  to  the 
source  duration  and  nonlinear  material  properties  will  not  be 
considered.  Their  inclusion  is  straightforward  and  does  not 
alter  the  conclusions  presented. 
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Consider  an  energy  deposition  profile  as  shown  in  Figure  D-2 
For  instantaneous  deposition,  the  initial  pressure  (determined 
by  the  Gruneisen  coefficient)  will  have  a  similar  form,  while 
the  initial  particle  velocity  is  zero  everywhere.  The  heated 
material  will  expand  both  at  the  free  surface  and  against  cooler 
material,  sending  a  compressive  stress  pulse  into  the  cold  region 
The  subsequent  deceleration  of  the  material  expanding  toward  the 
free  surface  produces  a  tensile  pulse,  if  no  spall  occurs. 

The  details  of  the  stress-wave  interactions  that  produce  the 
tensile  state  can  be  illustrated  by  considering  the  initial 
pressure  to  be  the  superposition  of  two  compressive  stress  pulses 
propagating  in  opposite  directions,  each  with  the  same  width  and 
half  the  amplitude  of  the  initial  pressure  distribution.  The 
tensile  pulse  is  thus  the  result  of  the  reflection  of  the  left¬ 
going  compressive  pulse  from  the  free  surface.  The  reflected 
pulse  can  be  treated  as  a  right-going  virtual  tensile  pulse. 

These  waves  are  shown  schematically  in  Figure  D-3. 

The  tension  develops  as  the  virtual  pulse  propagates  into 
the  material  (and  becomes  real) .  Since  the  tensile  pulse  follows 
behind  the  right-going  compressive  pulse,  it  is  superimposed  only 
with  the  left-going  compressive  pulse  (which  passes  out  of  the 
material  and  becomes  virtual) . 

First,  examine  the  case  where  the  spall  strength  of  the 
material  is  larger  than  the  peak  tension.  The  development  of 
the  stress  pulse  in  this  case  is  illustrated  sequentially  in 
Figure  D-4.  There  it  may  be  seen  that  the  peak  tension  at  each 
instant  occurs  at  the  location  of  the  leading  edge  of  the  tensile 
pulse  and  with  a  magnitude  determined  by  the  sum  of  the  maximum 
tensile  pulse  amplitude  and  the  magnitude  of  the  left-going 
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Figure  D-4b  Development  of  stress  waves  produced  by  sudden 
energy  deposition,  without  spall. 
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compressive  pulse  at  that  point.  Thus,  the  envelope  of  peak 
tension  is  formed  as  the  maximum  amplitude  of  the  tensile  pulse 
is  superimposed  successively  with  various  points  along  the  left¬ 
going  compressive  pulse.  Hence,  this  envelope  will  have  the 
same  shape  as  the  left-going  compressive  pulse,  except  that  the 
scale  of  the  abscissa  will  be  reduced  in  half  because  the  waves 
are  moving  toward  each  other  (Figure  D-5) .  Once  the  leading 
edge  of  the  tensile  pulse  has  reached  the  end  of  the  left-going 
compressive  pulse,  there  can  be  no  further  increase  in  the  peak 
value  of  the  tension.  This  occurs  at  tne  point  midway  between 
the  front  surface  and  the  end  of  the  region  of  energy  deposition, 
i.e. , half- range. 

Next,  consider  the  case  where  the  peak  tension  exceeds  the 
spall  strength.  For  instance,  if  spall  occurs  before  the  tensile 

stress  reaches  its  maximum  value,  as  shown  in  Figure  D-6 ,  a  new 

* 

free  surface  is  suddenly  created,  requiring  a  zero-stress 
boundary  condition.  At  this  point  the  situation  is  similar  to 
the  initial  condition.  Additional  tension  results  from  reflec¬ 
tion  of  the  left-going  compressive  pulse  remaining  in  the  un¬ 
spalled  material.  The  envelope  of  peak  tension  for  this  pulse 
may  be  constructed  as  in  the  first  case;  the  shape  of  the  curve 
is  unaltered  from  the  form  obtained  when  no  spall  was  allowed, 
except  that  it  is  translated  upwards  by  an  amount  equal  to  the 
spall  strength  (Figure  D-7) . 

If  the  magnitude  of  the  tensile  stress  pulse  is  more  than 
twice  the  spall  strength,  multiple  spalls  will  result.  The  en¬ 
velope  of  peak  tension  can  be  constructed  as  before,  and  the  last 


In  reality,  solid  spall  occurs  over  a  period  of  time  and  is 
more  complex  than  the  above  discussion  might  imply. 


157 


spall  occurs  when  the  remaining  left-going  compressive  pulse 
after  reflection  is  not  large  enough  to  exceed  the  spall  strength 
(Figure  D-8)  .  Thus,  even  if  a  zero  spall  strength  is  assumed, 

it 

no  front-surface  spall  occurs  beyond  half  range  since  there  is 
no  more  left-going  compressive  stress  pulse  remaining  to  reflect 
and  produce  tension. 

These  observations  apply  specifically  to  front-surface-spall 
phenomena.  In  making  comparisons  to  experimental  data,  it  should 
be  borne  in  mind  th-'t  the  right-going  compressive  pulse  may  re¬ 
flect  off  a  free  rear  surface  and  return  to  the  front  surface  as 
a  left-going  tensile  wave.  Thus,  unless  precautions  (such  as 
the  use  of  momentum  traps)  are  taken  to  prevent  this  complication, 
additional  material  loss  may  occur. 

Including  a  finite  source  time  reduces  the  amplitude  of  the 
peak  propagating  compressive  stress  and  peak  propagating  tensile 
stress  but  does  not  change  the  locations  of  the  maxima. 

At  this  point,  an  explanation  of  the  appearance  of  the 
craters  becomes  possible.  At  low  electron  energies  and  low 
fluence,  the  material  loss  is  due  to  either  solid  spall  or  a 
combination  of  liquid  spall  at  the  surface  and  solid  spall  deeper 
into  the  material.  In  all  such  cases,  the  crater  depth  is  less 
than  or  equal  to  the  half  range  of  the  energy  deposition.  Since 
a  substantial  gradient  in  temperature  is  obtained  between  the 
front  surface  and  the  half  range  point,  the  spall  strength  would 
increase  with  increasing  depth  into  the  material.  Hence,  it  is 

* 

In  performing  computations  using  finite-difference  material- 
response  codes,  care  must  be  taken  to  ensure  that  computational 
noise  does  not  cause  spall  indications  when  low  values  for  the 
threshold  are  used. 
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not  surprising  that  at  low  fluences,  the  spall  depths  are  less 
than  half  range,  since  the  material  at  the  half  range  depth  may 
be  able  to  sustain  considerable  tension. 

As  the  fluence  increases,  the  solid-partial  melt  interface 
moves  deeper  into  the  material.  As  it  passes  the  half  range 
point,  liquid  spall  occurs  only  up  to  the  half  range  and  melted 
material  is  left  on  the  sample. 

Further  increases  in  fluence  lead  to  larger  thicknesses 
of  melt  remaining.  The  radial  flow  patterns  observed  in  such 
cases  are  the  combination  of  radial  stress  relief  at  the  edge 
of  the  irradiated  area  and  small  variations  in  beam  intensity 
across  the  diameter  of  the  sample. 

When  the  higher-energy  beam  was  used,  melt  was  not  achieved 
to  depths  near  half  range.  Hence,  the  last  spall  was  a  solid 
spall,  leaving  a  clean  crater. 

The  above  discussion  provides  a  consistent  explanation  of 
the  observed  effect.  The  physical  phenomena  are  somewhat  differ¬ 
ent  when  vaporization  occurs.  If  extensive  vaporization  is 
produced  at  the  surface,  the  high-pressure  vapor  changes  the 
boundary  condition  at  the  liquid-vapor  interface.  As  a  result, 
the  tensile  pulse  may  be  suppressed,  and  little  or  no  liquid 
(or  solid)  spall  will  occur. 

The  situation  is  more  complex  when  a  substantial  portion  of 
the  irradiated  material  is  in  the  liquid-vapor  mixed  phase  region. 
In  such  cases,  it  appears  that  a  liquid  spall  might  occur,  while 
a  low-pressure  vapor  remains  in  front  of  the  spall  region,  in 
which  case  breakup  of  the  liquid  spall  could  allow  late-time 
communication  between  the  vapor  and  the  sample  surface. 
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Finally,  it  should  be  noted  that  the  total  impulse  is  com¬ 
paratively  insensitive  to  the  exact  quantity  of  mass  removed  by 
liquid  spall.  The  final  stages  of  liquid  spall  occur  in  material 
with  low  particle  velocity  so  that  the  impulse  associated  with 
this  mass  is  small. 
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ELECTRON  DEPOSITION  PROFILES  IN  ALUMINUM 
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Fiyure  Ii-3  Normalized  exposition  profile  in  aluminum. 


1/gram/cal/cm 


1/gram/cal 


1/gram/cal/cm 


